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ABSTRACT 

We present a new mass estimate of a well-studied gravitational lensing cluster, Abell 1689, from 
deep Chandra observations with a total exposure of 200 ks. Within r = 200 h^^kpc, the X-ray mass 
estimate is systematically lower than that of lensing by 30-50%. At r > 200 /i~^kpc, the mass density 
profiles from X-ray and weak lensing methods give consistent results. The most recent weak lensing 
work suggest a steeper profile than what is found from the X-ray analysis, while still in agreement 
with the mass at large radii. Fitting the total mass profile to a Navarro-Frenk- White model, we 
find M200 = (l-16to;27) X IO^^/i-^Mq with a concentration 



C200 



5.3^i['2, using non-parametric 



mass modeling. With parametric profile modeling we find M200 = (0.94^Qgg) x lO^^/i ^M0 and 
C200 = 6.6lg'4. This is much lower compared to masses deduced from the combined strong and 
weak lensing analysis. Previous studies have suggested that cooler small-scale structures can bias 
X-ray temperature measurements or that the northern part of the cluster is disturbed. We find these 
scenarios unlikely to resolve the central mass discrepancy since the former requires 70-90% of the space 
to be occupied by these cool structures and excluding the northern substructure does not significantly 
affect the total mass profiles. A more plausible explanation is a projection effect. Assuming that the 
gas temperature and density profiles have a prolate symmetry, we can bring the X-ray mass estimate 
into a closer agreement with that of lensing. We also find that the previously reported high hard- 
band to broad-band temperature ratio in A1689, and many other clusters observed with Chandra, 
may be resulting from the instrumental absorption that decreases 10-15% of the effective area at ~ 
1.75 keV. Caution must be taken when analyzing multiple spectral components under this calibration 
uncertainty. 

Subject headings: galaxies: clusters: individual: Abell 1689 — X-rays: galaxies: clusters 



1. INTRODUCTION 

Abell 1689 is a massive galaxy cluster with the largest 
know n Einstein radius to date, 9e = 45" for Zs = 1 



affect lensing mass estimates. 
A recent joint Chandra, 



HST/ACS, an d Sub 



(2008a 



1^ 



Tyson et al."1990l 



Miralda-Escude & Babul 1995 



Broadhurst et al. 2005a Tbt, located at a moderately low 



0.187 (I'rye et al.|[2007[ ). It has a regular 



redshift of z 

X-ray morphology, indicating that the cluster is likely in 
hydrostatic equilibrium, but the mass derived from the 
X-ray measurement is often a factor of 2 or more lower 
than that from gravitational len sing at most radii. Us- 
i ng X MM-Newton observations, [Andersson fc Madejski| 
(2004 A04 hereafter) find an asymmetric temperature 
distribution and a high redshift structure in A1689, pro- 
vi ding evidence for a n ongoing merger in this cluster. 



aru/Suprime cam analysis by Lemze et al 
L08 hereafter) suggested that the temperature of A168!) 
could be as high as T = 18 keV at 100 /i~^kpc, almost 
twice as large as the observed value at that radius. The 
derived 3D temperature profile was based on the X-ray 
surface brightness, the lensing shear, and the assumption 
of hydrostatic equilibrium. From the disagreement be- 
twee n the observed X-ray temperature and the deduced 
one, |L08| concluded that denser, colder, and more lumi- 
nous small-scale structures could bias the X-ray temper- 
ature. 

In another stud y of 192 clusters o f galaxies from the 



Saha et al.| ( 2007 1 confirm the existance of substruc 
tures, using diilerent sets of lensed imag es. This is also 



Chandra archive, Cavagnolo et al. ("2008) find a very 



seen in other lens i ng work (e.g., fe roadhur st et al.|2005a 
Diego et a l. 2005 : Z gkser et al7|[2P06[ .Halkoia ei alM im 
Limousin et al. 2007|. Though these clumps are cli 
ly ' 



high hard-band {2/[l+z)-7 keV) to broad-band (0.7-7 

compared 



1.36 



-1-0.14 
-0.12 



clumps are clearly 

identified, they only^ contribute ~ 7% of the total mass 
within 250 ft.~^kpc and are likely to be line-of-sight fila- 
ments rather than distinct merging groups. Furthermore, 
Lokas et al. ( 2006 ) used the redshift distribution of galax- 
ies to conclude that A1689 is probably surrounded by a 
few structures superposed along the line of sight that 
do not interact with the cluster dynamically, but would 



keV) temperature ratio for A1689, 
to 1.16 ± 0.10 for the whole sample. They also find that 
merging clusters ten d to have a higher tem perature ra- 



tio, as predicted by Mathiesen & Evrard (j2001J where 



this high ratio is attributed to accreting cool subclus 
ters lowering the broad-band temperature by contribut- 
ing large amounts of line emission in the soft band. The 
hard-band temperature, however, s hould be unaltered by 
this em ission. The simulations of IMathiesen k, EvrardI 
(2001 1 show an increase of temperature ratios of ~ 20% 
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in general, which is clo se to the average of the sample of 
Cavagnolo et al.lpOOSl), 16%. 



A recent study, us ing the latest Chandra data (Riemer 



S0rensen et al. 2009 1 claim that the cluster harbors a cool 



core and thus is relaxed based on a hardness-ratio map 



2 



analysis. They further calculate a mass profile from the 
X-ray data and conclude that the X-ray and lensing mea- 
surements are in good agreement when the substructure 
to the NE is excluded. 

In this work, we examine the possibility of an extra 
spectral component in the X-ray data and derive an im- 
proved gravitational mass profile, including a recent 150 
ks Chandra observation. S|2] describes the details of data 
reduction and examines the possibility of an uncorrected 
absorption edge in the data. In f|3] wc explore the physi- 
cal properties of the potential cool substructures under a 
two-temperature (2T) model and examine if they can be 
used to explain the high hard-band to broad-band tem- 
perature ratio. I n ^ assuming that the temperature 
profile derived by [L08' is real, we investigate what this 
implies for the required additional cool component. In 
^ we derive the mass profile under both one and two 
temperature-phase assumptions, using both parametric 
and non-parametric methods. Finally, we discuss our re- 
sults in fp] and summarize in ^ 

Throughout this paper, we assume Hq = 100 h^^ km 
s~^ Mpc~^, ^Im = 0.3, and fi^ — 0.7, which give s 1" = 



2.19 h ^kpc at the cluster redshift of 0.187 (Frye et al. 
2007|. Abund ances are relative to the ph otospheric solar 



All errors are 



abundances of Anders & Grevesse ( 1989 1 
1 a unless otherwise stated. 

2. DATA REDUCTION 
2.1. Chandra 

Chandra data were processed through CIAO 4.0.1 with 
CALDB 3.4.3. Since all of the observations had gone 
through Repro III in the archive, reprocessing data was 
not needed. Updated charge-transfer inefficiency and 
time-dependent gain corrections had already been ap- 
plied. For data taken in VFAINT telemetry mode, addi- 
tional screening to reject particle background was used. 
Events with bad CCD columns and bad grades were re- 
moved. Lightcurves were extracted from four I-chips 
with cluster core and point sources masked in the 0.3- 
12 keV band and filtered by Ic-clean which used 3 a 
clipping and a cut at 20% above the mean. Finally, 
make-readout^hg were used to generated Out-of-Timc 
event file. These events were multiplied by 1.3% and sub- 
tracted from the images or the spectra to correct read- 
out artifacts. For spectral anlysis, emission- weighted re- 
sponse matrices and effective area files were constructed 
for each spectral region by mkacisrmf and mkwarf. 

2.2. Background Subtraction and Modeling 

Blank-field data sets'^ were used to estimate the back- 
ground level. After reprojecting the blank-sky data 
sets onto the cluster's sky position, the background was 
scaled by the count rate ratio between the data and the 
blank- field background in the 9.5-12 keV band to account 
for the variation of particle induced background. Below 
1 keV, the spatial varying galactic ISM emission (Marke^ 



vitch et al.| 2003 1 could cause a mismatch between the 
real background and the blank-field data. By analysing 
the spectra in the same field but sufficiently far from the 
cluster, tailoring this soft component can be made using 

3 |http : //asc . harvard ■ edu/cal/Acls/Cal_prods/bkgrnd/| 

|acisbg| 
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Fig. 1.— The 0.6-9.5 keV Chandra spectrum of A1689 from the 
central 3' region. The upper panel shows the data, plotted against 
an absorbed VAPEC model (solid line) with each element's abun- 
dance and absorption column density as free parameters. The lower 
panel shows residuals. 

an una bsorbed T ^ 0.2 keV , solar abundance thermal 
model dVikhlinin et al.||2005[ ). 

The current available blank-sky data were created from 
observations before 2005. As the solar cycle gradually 
reaches its minimum, the particle induced background 
increases. Therefore, newer observations need a much 
higher background normalization with a factor of 1.2- 
1.3. This leads an overestimate of the background in 
the soft band because other components like cosmic X- 
ray background (CXB) does not change as the particle 
induced background does. To correct the over-subtracted 
CXB and halo emission, an absorbed power law with 
photon index fixed at 1.4 ( De Luca fc Molendi|2004 ) plus 
an unabsorbed thermal model was used to tit the blank- 
field background subtracted spectrum taken at r > 13'. 

The background normalization factors used for each 
observation are hsted in Table[T] 



2.3. XMM-Newton 

The data from two MOS detectors were processed with 
the XMMSAS 6.1.0 tool, emchain. Backgr ound flares 



were remov ed by a double- filtering method (Nevalainen 
et al.||2005l) from E > 10 keV and 1-5 keV light curves 



Only events with pixel PATTERNS 0-12 were selected. 
Since XMM data were only used to crosscheck the re- 
sult of the multi-component analysis of Chandra spec- 
tra, extracted from the central region where background 
modeling is relatively unimportant, we used the simpler 
local background, taken from 6'-8'. Spectral response 
files were created by rmfgen and arfgen. We did not 
include PN data because the measured mean redshift, 
0.169 ± 0.001, was not consistent with those from XMM 
MOS or Chandra data. This could i ndica te a possible 
gain offset for PN detector, although A04 did not find 
any evidence for that. 

2.4. Systmatic Uncertainies 

|L08| pointed out some issues about previous Chandra 
observations (ObsID 540, 1663, and 5004). The column 
density from Chandra data is much lo wer than the Galac- 
tic value, 1.8 x 10^° cm"^ (Dickey fc Lockman||r990| ), 
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TABLE 1 
Chandra Observation Log 



ObsID 


Data 


Obs. Date 


Exposure 


Back J 


;round Normalization 




Mode 




(ks) 


10 


11 


12 


13 


540 


FAINT 


2000-04-15 


10.3 


1.06 


1.09 


1.03 


1.11 


1663 


FAINT 


2001-01-07 


10.7 


1.00 


0.98 


0.99 


1.04 


5004 


VFAINT 


2004-02-28 


19.9 


0.93 


0.89 


0.89 


0.94 


6930 


VFAINT 


2006-03-06 


75.9 


1.21 


1.18 


1.18 


1.28 


7289 


VFAINT 


2006-03-09 


74.6 


1.19 


1.20 


1.19 


1.27 



10 

-10 
-20 



10- 



-10 
10 



-10- 




-ObsID 6930+7289 
-ObsID 5004 
ObsID 540+1663 




Energy (keV) 

Fig. 2. — Fit residuals, showing each channel's contribution to the 
total x^- Top: an absorbed VAPEC model fit to central 3' spec- 
trum. Middle: same as the above, ignoring data in 1.7-2 keV. 
Bottom: adding an absorption edge with E-tii,.^^}^ =1.77 keV and 
r=0.12. 



which is also supp orted by the ROSAT data (Andersson 
fc: Madejski|2004 ) . The temperature difference can be as 
high as 1.'6 keV depending on the choice of column den- 
sity. In the high energy band, the data is systematically 
higher than the model prediction. With two long Chan- 
dra observations, ObsID 6930 and 7289, we clearly see 
an unusual feature in the datasets which may give clues 
to problems r nentioned before. Fig.[l]shows an absorbed 



APEC model ( [Smith et al. 12001 1 ) fit to the central 3' spec- 
trum. The prominent residual at '--^ 1.75 keV is present in 
all of our observations and appeared as the biggest con- 
tributor to the total (See Fig. 2]). This residual can 
not be eliminated by adjusting individual abundances in 
the cluster or in the absorbing column (the cluster is at 
high galactic latitude). Because the residual around 1.75 
keV is an order of magnitude larger than the background, 
it is not likely related to the background subtraction. In 
addition to this absorption, the residuals are systemat- 
ically rising with the energy from negative to positive 
values. This trend is not changed when fitting the spec- 
trum with data between 1.7-2.0 keV excluded (Fig. [2]). 
We found that multiplying a XSPEC Edge model can 
correct the residual at ~ 1.75 keV, remove steadily rising 
residuals with the energy, and make the column density 
agree with the Galactic value. 



TABLE 2 

Absorption edge parameters 



Model 


fit range 


Ethresh 


T 






(keV) 





IT 

IT 

IT 

2T^ 

ITH 

IT 

IT 



PLt- 



2.5'x2.5' 1.75 

■r <3' 1.74; 

r <3', ignore 1.75-1.85 keV 1.76 

■r <3', ignore 1.75-1.85 keV 1.76 

r <3', ignore 1.75-1.85 keV 1.76 

0.2'< r <3' 1.74; 
0.2'< r <3', ignore 1.75-1.85 keV 1.90 



+0.01 
0.01 
+0.01 
0.01 
+0.03 
0.02 
+0.06 
0.02 
+0.03 
0.02 
+0.01 
0.01 
+0.04 
0.06 



0.15 
0.14 
0.13 
0.12 
0.11 
0.14 
0.12 



+0.01 
0.01 
+0.01 
0.01 
+0.01 
0.01 
+0.02 
0.01 
+0.01 
0.01 
+0.01 
0.01 
+0.01 
0.01 



S.Otg J keV, T2 = keV. 



: 9.3^ 



keV, r : 



-0.7 



+0.3 
0.4- 



Since the spectrum was extracted from a very large re- 
gion, we averaged the position-dependent response ma- 
trices and effective area functions by the number of 
counts at each location. It is possible that the absorp- 
tion feature is caused by improper weighting of those 
response files, or that this peculiarity only exists at cer- 
tain regions. To dispel those doubts, we separated the 
central 2.5' x 2.5' area into 12 square regions and simul- 
taneously fit these spectra with one spectral model (We 
only used data from ObsID 6930 and 7289 to simplify 
the fitting procedure). All parameters, except for the 
normalization, were tied together. The residuals from 
the single temperature fit are shown in Fig. [3] Although 
the fit is now acceptable with a /dof = 3448.8/3387, 
the residuals still show the same systematic trend as 
seen in the composite spectrum, and the column den- 
sity, (l.Ot^j) X 10^0 cm-2 (90% confidence level), ap- 
pears low. When adding an absorption edge to the single 
temperature model, the derived parameters of this edge, 
Ethresh = l-75l°:oi keV and r = 0.15l:o;°i, are consis- 
tent with results from the integrated spectrum. In fact, 
Ethresh and r do not strongly depend on how we model 
the cluster spectrum. We list fitted values of Ethresh 
and r from different cluster models and spectral extrac- 
tion regions in Table[2j Similar values are also found in 
other Chandra datasets (see the Appendix) . 

The low column density can be explained by the ab- 
sorption at ~ 1.75 keV. This absorption has to be some- 
thing like an edge, which affects the wide energy range of 
the spectrum, rather than an absorption line that influ- 
ences only a limited energy range. The consequence to 
the fit resulting from this absorption is that the overall 
spectral normalization decreases. Since less soft photons 
are produced from the model, the heavy absorption by 
the foreground is no longer needed. We quantified the 
above statement in Table|3] in which we simulated 100 
single temperature spectra with an absorption edge at 



4 




Energy (keV) 

Fig. 3. — Residuals from an absorbed APEC model fit to 12 
spectra extracted from central 2.5' X 2.5' region of ObsID 6930 and 
7289. We simultaneously fit these spectra and tied all parameters, 
except the normalization, together. 

Ethresh = 1-74 keV and t = 0.14, fit with a single tem- 
perature (IT) model, and compared with the observa- 
tions. The spectral normalization is increasing with ris- 
ing column density as we exclude more data around 1.75 
keV. Meanwhile, the cluster temperature and abundance 
are slightly decreasing. The changes of those parame- 
ters from different bandpass used in the fitting match 
perfectly to what are seen in the real data. 

The CCD calibration around the Si-edge for ACIS-I 
detectors is a known issue (N. Schulz, private communi- 
cation). However, it is unknown whether a correction like 
an edge model is needed, or if we should simply ignore 
the data around the Si edge. If the former is true, re- 
sults from the multiple-component analysis or the hard- 
band/broad-band temperature measurement without ap- 
plying this edge model beforehand are very questionable. 
As seen in Table|3] including the edge model can make 
the hard-band temperature 30% hotter than the that of 
the broad-band. This temperature ratio depends on the 
cluster temperature and the quality of the data. On the 
other hand, if the latter is true, the spectrum implies that 
an additional component which is much harder than 10 
keV emission is definitely required. Though, the fit is 
not as good as that with an edge model. From the fact 
that the absorption depth is sufficiently far from the zero, 
even if we exclude data around 1.75 keV (Table[2]), the 
intensity jump around this energy indeed exists. 

A1689 is a very hot cluster that unfortunately will be 
seriously affected by the calibration uncertainty around 
the Si edge if that can be modeled by something similar 
to an edge model. Lacking the knowledge that correctly 
treats the systematic residuals seen in the data, we pro- 
vide both models, applying an absorption edge or simply 
ignoring the data around 1.75 keV, as our best guess to 
the thermal state of this cluster. 

2.5. Spectral fitting 



Spectra were fit with XSPEC 12.3.1 package (Arnaud 



We adopted x statistic and grouped the spec- 

tra to have a minimum of 25 counts per bin. However, 
when fitting background dominated spectr a, \^ statis 



tic is proven to give biased temperature (Leccardi & 
Molendi|2007l ). Another choice available in XSPEC is us- 
ing Cash statistic with modeled, rather than subtracted 
background. Since modeling the background needs many 
components: cosmic ray induced background (broken 
power-laws plus several Gaussians), particle background 
(broken power-laws), cosmic X-ray background (power- 
law) , galactic emission (thermal) , etc, the whole spectral 



model will be very complicated for analysis like depro- 
jection, which simultaneously fits all of the spectra ex- 
tracted at different radii. We decided to use statistic 
but with a different grouping method to bypass the dif- 
ficulty in background modeling. 

As shown in Tablej4] we simulated 500 Chandra spectra 
with Nh = 1.8x lO^^cm-^ Z = 0.2 solar, and T = 9, 7, 5 
keV with Norm^ = 1.51 x 10'*,1.52 x 10^,1.56 x lO"*, 
respectively. The spectral normalization was chosen to 
match the observed fiux at r = 6.5'-8.8', where the back- 
ground is ^ 90% of the source in 0.9-7.0 keV band. Spec- 
tra were generated based on the response files of ObsID 
6930 and 7289 with a total exposure time of 150 ks. 
When data are binned to have a minimum of 25 total 
counts (background included) per channel, a 9 keV gas 
will be measured to be 5 keV. Raising the threshold can 
lessen this bias. However, even with 400 counts per bin, 
which greatly reduces the spectral resolution by a factor 
of 10, the temperature is still being underestimated by ~ 
1 keV for T — 9 keV gas. We found out that binning data 
to have at least 2 counts above the background can re- 
cover the true temperature, though this minimum have 
to be adjusted according to the background contribu- 
tion. Spectra at large radii were binned by this grouping 
scheme. 

3. SPECTRAL ANALYSIS 
3.1. Single temperature model 

We first determined the general properties of the clus- 
ter, using the spectrum extracted from the central 3' (395 
h~^kpc) region and fitting it with a single temperature 
VAPEC model. The Ne, Mg, Si, S, Ar, Ca, Fe, Ni abun- 
dances and the redshift were free to vary. The column 
density was fixed at the Galactic value. The best-fit pa- 
rameters are listed in Table[5l For Chandra data we fixed 
t he Si abundance at 0.4 Zq, since the residual at 1.75 keV 
(p^Il was close to Si XIV Ka line (2.01 keV, rest-frame). 

i'ablejS] shows that a single-temperature model is mod- 
erately adequate for XMM-Newton MOS data but not 
for Chandra when the absorption edge is not modeled. 
In addition to this difference, the Chandra temperature 
is --1 keV higher than that of XMM-Newton. This 
temperature disagreement is likely related to the c ross 
calibration problems, as noted in other studies 



Ko- 



tov fc Vikhlmm '2005^ 'Vikhlinin ct al. 'MM, [Snowmen 
et al. 12008 L. David^), but it could also be caused by in- 



correct cluster modeling. A two-temperature model will 



be investigated in ^3.2 

Chandra data also have a higher Fe and a much higher 
Ni abundance, resulting an unusually high Ni/Fc ratio 
of 7.5 ± 1.3 and 9.0 ±1.2 Nig/Feo with and without 
an absorption edge correction, respectively, in contrast 
to the XMM-Newton value of 1.8 ± 1.6 Nio/FeQ. Our 
XMM-Newton MPS result is in agreement with that of 
de Plaa et al] ( |2007[ ) , 0.9 ± 1.5 Nio/Feo, obtained from 
MOS and PIN spectra from the r < 1.3' region with a dif- 
ferential emission measure MEKAL-based model, wdem 



( [Kaastra et al. 2004). Such a high Ni/Fe ratio greatly 



Spectral normalization, Norm 



10" 



47r{{l + z)DA)2 

where rie and nn are in cm~^, V in cm^ , and in cm. 

^ http : // cx c ■ harvard ■ edu/cc«/proceedings/07-proc/| 

present at ions/david 
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TABLE 3 

Summary of r < 3' spectral fits and IT simulations 









Chandra observations 






Simulated 200 ks spectra'' 


fit range 


T 


Z 


Nh Norm 


xVdof 


T 


Z Nh Norm 




(keV) 


(Zq) 


(102"cm-2) (10-2) 




(keV) 


(Zq) (1020cm-2) (10-2) 



0.6-9.5 keV*" 10.5 

0.6-9.5 keV 10.7; 

ignore 1.75-1.85 keV 10 

ignore 1.7-2.0 keV 10.5 

ignore 1.7-2.5 keV 10.4 

ignore 1.7-3.0 keV 10.4 

1.7-7.0 keV 13.1 

0.7-7.0 keV 10.1 



+0.1 
0.1 
+0.4 
0.1 
+0.1 
-0.1 
+0.1 
0.1 
+0.1 
0.1 
+0.1 
0.1 
+0.3 
0.3 
+0.1 
0.1 



0.36 
0.41 
0.40 
0.39 
0.38 
0.37 
0.44 
0.40 



+0.02 
0.02 
+0.02 
0.02 
+0.02 
0.02 
+0.02 
0.02 
+0.02 
0.02 
+0.02 
0.02 
+0.02 
0.02 
+0.02 
0.02 



1.9: 

0.7l 
0.9"* 



3+0.2 
-0.2 
7+0.2 
-0.2 
3+0.2 
-0.2 

1 l+"-2 
^■^-0.2 
„+0.2 
■-'-0.2 
=+0.2 
'-0.2 

.8 



1 920+°°°'' 
1 SA1 +0-007 

1.0'il_g QQy 

-, or;7+0.007 

1 SKS+0-008 
1-»D»_0.008 
1 SQ9+0-008 

1 904+0 009 

1.3U1_o QQg 

1 o-i 7+0.008 
^■"^ ' -0.008 
1 Sfi7+0-006 
1-»D'_0.006 



1591/1390 
1867/1390 
1743/1366 
1645/1327 
1466/1222 
1317/1120 
1100/960 
1655/1164 



10.4l°;l 
10.6l°:l 
10.5l°:l 

io.3i°:l 

10.2l°;l 

io.2l°;i 

-19 7+0.3 
-^-^•'-0.3 

io.oi°:l 



0.36 
0.40 
0.39 
0.38 
0.37 
0.36 
0.44 
0.39 



02 
02 
02 
02 
02 
02 
02 
02 
02 
.01 
.02 
01 
02 
.02 
.02 
02 



1 




1.2 

1. 

1. 

1. 

1. 



+0.1 
-0.2 
+0.1 
-0.2 
+0.1 
-0.1 
+0.2 
0.1 
+0.2 
-0.2 
+0.2 
-0.2 



1 899+0-008 

1 090+0-008 
i.»^0_o 007 

-I OQQ+0-008 

l.i500_o 006 
1 sc:9+0.008 
±.OOZ_o 007 
1 Q77+0-007 
' -0.008 
1 SSfi+0-008 
i.»»6_o 009 
1 700+0-008 
l.l»u_o 008 
1 849+0-005 
l.O'i»_o 006 



Note. — Errors arc Icr or 68% CL for 100 simulations. 
' The absorption edge is at Ethresh — 1.74 keV with r — 0.14. 



to Sx[0.6-9.5fceV] = 2.4 X 10 
b 



Wo used T -. 

-2 



erg cm-^s-^), Nh = 1.8 X 10^" cm 
Multiplied by an absorption edge at Ethresh — 1.74 keV with r — 0.14 



TABLE 4 

Summary of IT simulations of 150 ks Chandra 
SPECTRA at r = 6.5'-8.8'. 



Min counts 
(1) 


To (keV) 
(2) 


T„ed (keV) 
(3) 


do f mean 

(4) 


25 (tot) 


9 
7 
5 


4.9t2-i 
3.4t°:? 


396 
396 
395 


100 (tot) 


9 
7 

5 


5.9t?:i 


142 
142 
141 


400 (tot) 


9 
7 

5 


7Q+7.2 

6.6t^:2 
4.8tl:«3 


40 
40 
40 


2 (net) 


9 
7 
5 


8.6+^-^7 
7 2+''-^ 

'•^-2.2 
5 9+1-9 


170 
166 
155 



Note. — (1) minimum (total or net) counts per 
channel, (2) input temperature, (3) median temper- 
ature, (4) median dof. Errors are 68% CL for 500 
simulations. 



excee ds the yield of typical SN la models ( Iwamoto et al. 
19991, which range from 1.4-4.8 Nio/Fe©. Since there is 
a known temperature discrepancy between Chandra and 
XMM-Newton that would affect elemental abundance de- 
terminations, direct Fe and Ni line measurements will be 
conducted in 3 3. 4| 

The S abundance, determined mostly by the S xvi Ka 
line at 2.62 keV (rest-frame), should be accurately mea- 
sured for XMM-iVe wton EPIC since it suffers little sys- 
tematic uncertainty (Werner et al. 2008). However, with- 
out an absorption edge correction, there is basically no 
S detection for Chandra data, which strongly contradicts 
the XMM-Newton result. This shows the great impact 
of the absorption at 1.75 keV. When including an edge 
model into the fit, we have consistent S abundances for 
both instruments. 

3.2. Two temperature model 

To get some clues to the nature of the claimed cool 
substructures in A1689, a simple two temperature model 
was fit to the spectrum extracted from the r < 3' (395 



10.5 keV, Z = 0.36 Zq, z = 0.187, Norm = 1.9 X 10"^ (correspondii 



/i~^kpc) region where the quality of the data was high 
enough to test it. We used two absorbed VAPEC mod- 
els, with variable normalization but linked metallicities 
between the two phases. The column density was fixed at 
the Galactic value. To reduce the uncertainty on measur- 
ing metallicities, we tied the abundances of a-elements 
(O, Ne, Mg, Si, S, Ar, and Ca) together and fixed the 
remaining abundances at the solar value, except for Fe 
and Ni. Since the hotter phase temperature, Thot, was 
harder to constrain, it was frozen at a certain value above 
the best-fit single temperature fit. Tit- We changed this 
increment from 0.5 to 50 keV to explore the whole pa- 
rameter space. 

Fig.|4] shows the temperature of the cooler gas, 
Tcooi, and the fractional contribution of the cooler gas, 
EMcooi/EMtotai, as a function of Thot- As Thot increases, 
Tcooi and EMcooi/EMtotaiincrease as well. Tcooi eventu- 
ally becomes Tit once Thot is greater than 20 keV and 
very little gas is left in the hot phase, which is also sup- 
ported by the XMM-Newton data. For Thot « 18 keV, 
there has to be 30%, 60% of the cool gas at the tempera- 
ture of 5, 8 keV inferred from Chandra and XMM-Newton 
data, respectively. Chandra absorption corrected data 
show similar results as XMM-Newton data do at this tem- 
perature. Although there is some inconsistency between 
Chandra and XMM-Newton data, both indicate that the 
cool component, if it indeed exists, is not cool at all. 
T — 5 keV is the typical temperature of a medium sized 
cluster with a mass of Afsoo = 2.9 X lO^-^/t-^Me (Vikhlini^ 



et al.|20 06). 

To quantify how significant the detection of this ex- 
tra component was, we conducted an T-test from the 
fits of IT (the null model) and 2T models. However, 
because the 2T model reduces to IT when the normal- 
ization of one of the two components hits the parameter 
space bound ary (ie, zero), the assum ption of T-test is not 
satisfied (see Protassov et al.|2002 ). Therefore, we simu- 
lated 1000 11' Chandra spectra and performed the same 
procedure to derive the F-test probability, Pp, based on 
the F distribution. Fig. [4] shows the distribution of Pp 
from simulated data at the 68, 90, 95, and 99 percentile 
overplotted with Pp from Chandra and XMM-Newton 
data. We plot Pp in Fig. |4] rather than the F statistic, 
since Pp is a scaler that does not depend on the degrees 
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TABLE 5 
Best-fit VAPEC parameters 





T (keV) 


z 


Nc 


Mg 


Si 


s 


Ar 


Ca 


Fc 


Ni 


xVdof 


Chandra ^ 
Chandra ^ 
XMM MOS 


10.2«;1 
10.5«:1 
9.41°;} 


U.iSb^o.ooi 

U.i«ci_o 001 


U.OM_o 29 

< 0.18 

'^■'^0-0.35 


1 Q5+0.30 

1-1D_0.42 
l.OO_0.40 


0-4/ 
0.4 J 

0.501°;?° 


< 0.05 
0.66l«o-jl 

U.D'l_o 23 


< 0.19 

< 0.38 

< 0.57 


< 0.12 

< 0.09 

-I qo+0.69 


Q 07+0.01 
-0.02 

'^■'-'^-0.02 


■3.DU_o 47 
+0.48 
'-0.45 
U.O(_o 47 


1664/1360 
1545/1382 
982/840 



Note. — Al, O fixed at 0.4 solar and He, C, N at 1 solar. For the elements whose abundances reach the lower bound, zero, only the upper limits are 
shown. Errors arc Icr. 

^ without an absorption edge and ignoring data at 1.75-1.85 kcV. 



with an absorption edge. Edge parameters are determined from the 
T 1 1 1 



data with Ethrcsh = i-^ito m IfeV and r = 0.13t° qI- 
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Fig. 4. — The temperature of the cooler gas Tcooii the emission 
measure ratio EMcool/EMtotali a^nd the F-test probability Pp are 
plotted as a function of Tjjot. The shaded region represents 68%, 
90%, 95%, and 99% CL from 1000 simulated T = 10.1 keV Chandra 
spectra. The Pp from Chandra data without the absorption edge 
corrected (circles) is multiplied by 10^. 

of freedom of the fits and is ideal to compare observations 
that have different data bins. For Thot < 20 keV, both 
the edge-corrected Chandra data and the XMM-Newton 
data are within the 95 percentile of the simulated IT 
model and we conclude that a 2T model is possible but 
not necessary to describe the data. 

3.3. Hard-band, broad-band temperature 

In addition to multiple-component modeling, measur- 
ing the temperature in different band-pass is another 
way to dem onstrate the presence of multiple components. 
Cavagnolo et al. ( 20081) reported a very high hard-band to 
broad-band temperature ratio for A1689, 1.36^q ^[ji from 
analysis of 40 ks of Chandra data, suggesting that this 
could relate t o ongoing or recent me rgers. Following the 
convention in Cavagnolo et al. ( 2008 1, we fit the spectrum 
in the 0.7-7.0 keV (broad) and 2.0/(l-Fz)-7.0 keV (hard) 
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Fig. 5. — The hard-band to broad-band temperature ratio 
T2. 0-7. o/To. 7-7.0 of simulated Chandra 2T spectra (circles) plot- 
ted against Ti^of The shaded regions show the observed temper- 
ature ratios from Chandra and XMM-Newton MOS data. Also 
shown is the tempe rature ratio from 40 ks Chandra data by |Cav-] 
lagnolo et"ar] ( |2008^ . 



band with a single-temperature model. In contrast to 
C08, we do not use the r < i?2500 region with the core 
excised, but simply take the spectrum from the whole 
central 3' (395 /i^^kpc) region. The hard-band to broad- 
band temperature ratio from Chandra data, 1.29 ± 0.03, 
strongly disagrees with that of XMM MOS, 1.07 ± 0.03. 
This result is anticipated since an a bsorp tion edge fea- 
ture found in the Chandra spectrum (^ 2.4 1 is close to the 
cut-off of the hard band. After correcting lor this absorp- 
tion, the temperature ratio is in the range of 1 to 1.08 for 
an absorption depth of r = 0.14 — 0.10. As a consistency 
check, we simulated spectra accordi ng to the best-fit 2T 
models (from Chandra data) from S3.2 to see whether 
these models can explain such a high temperature ratio. 
Results are plotted in Fig.|5] None of the 2T models can 
reproduce the observed ratio of the uncorrected Chandra 
data. Thus we conclude that there is no evidence from 
this ratio of the presence o f multiple component s or me rg- 
ing activity. Furthermore, Leccardi & Molendi (20081 do 
not find any discrepancy between the hard band (2-10 
keV) and broad band (0.7-10 keV) temperature profiles, 
except for r < 0.05 rigo, for a sample of ~ 50 hot, inter- 
mediate redshift clusters based on XMM-Newton obser- 
vations. The high hard-band to broad-band temperature 
ratio seen in A1689, as well as in many other c lusters ob- 
served with Chandra (Cavagnolo et al. 20081, might be 
due to the aforementioned calibration uncertainty. 



3.4. Emission line diagnostics 

When fitting the whole spectrum, the temperature is 
mainly determined by the continuum due to the low 
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TABLE 6 



Line 


Energy 


Centroid ^ 


Width =>• 




(keV) 


(keV) 


(eV) 


Fe XXV Ka 


6.636, 6.668, 6.682, 6.700 


6.686 


23 


Fe XXV K/3 


7.881 


7.877 


19 


Fe 8.3 keV'= 


8.246, 8.252, 8.293, 8.486 


8.282 


68 


Fe XXVI Ka 


6.952, 6.973 


6.964 


14 


Fe 8.7 keV 


8.698, 8.701, 8.907, 8.909 


8.764 


97 


Ni XXVII Ka 


7.765, 7.805 


7.793 


19 


Ni XXVIII Ka 


8.074, 8.101 


8.090 


16 



TABLE 7 

Summary of line analysis 



Emissivity-weighted center and one standard deviation. The hne 
emissivity is calculated at T = 10 keV from Chandra ATOMDB 
1.3.1. 

^ including Fe xxvi K/3, XXV K7, and XXV Ki5. 
including Fe XXVI KS and xxvi K7. 

amount of line emission at the temperature of A1689. 
In order to extract the emission hne information, which 
can provide an additional temperature diagnostic, we 
fit the 4.5-9.5 keV spectrum with an absorbed ther- 
mal bremsstrahlung model plus Gaussians. There are 
42 lines whose emissivity is greater than 10~^^ pho- 
tons cm~'^s~^ at kT = 10 keV from ions of Fe xxv, 
Fe XXVI, Ni xxvii, Ni xxviii, according to Chandra 
ATOMDB 1.3.1. Considering the CCD energy resolu- 
tion, we grouped those lines into seven Gaussians and 
used the emissivity-weighted centroid and one standard 
deviation as the line center and width, respectively. The 
Ni XXVII Ka hne is ~ 80 eV away from the Fe xxv K/3 
line, not separable under CCD resolution unless we have 
extremely good data quality. Since we obtained an un- 
usually high Ni/Fe ratio of ~ 9 Nio/Fe© from a VAPEC 
model fit to the whole spectrum (j |3.1[ ), it is worth to in- 
vestigate this in detail. We therefore modeled Ni xxvii 
Ka and Fe xxv K/3 lines individually. Fig. ^ shows the 
spectrum and the best-fit model. The modeled lines are 
Usted in Tabled 

Strictly speaking, using fixed values of line centroids 
and widths is not correct because those quantities change 
with temperature. In addition, we approximated the line 
complex as a Gaussian whose line centroid and width 
calculated from the model may not be the same after 
being convolved with the instrument response. To prop- 
erly compare our fit results with the theory, we simulated 
spectra and fit them the same way we fit the real data. 
Fig. [7] shows the observed line ratios and results from 
simulated VAPEC spectra with 9 Nio/Fe©. 100 spectra 
were produced at each temperature and the flux was kept 
at the same level as that of the data. From the good 
match of fitted results from simulations to the direct 
model prediction, we confirmed that the fitting is accu- 
rate enough to measure the line flux, though only Fe xxv 
Ka and Fe xxvi Ka lines are precise enough for temper- 
ature determination. Table[7] shows the temperature and 
abundances, inferred from a single-temperature APEC 
model. The iron line temperature is in very good agree- 
ment with the continuum temperature for both Chandra 
and XMM-Newton data. All the Chandra and XMM- 
Newton observed line fluxes, except Fe xxv K/3, are con- 
sistent with each other (after an overall 9% adjustment 
to the flux). Using Fe xxv-t-xxvi Ka and Ni xxvii Ka 
line flux, we obtain accordant Fe and Ni abundances 
from both instruments. The larger Fe and Ni abun- 







Chandra 


XMM MOS 


Continuum 








T 


(kcV) 






Emission lines 








^ a 

Ni/Fe ''t 
Ni/Fe 


(keV) 

(Nio/Feo) 
(NiQ/Fee) 


9 6+° '^ 

8 4+3 '^ 
5 ci+3.2 


10 i+°-'^ 

1 4+1-8 

Q 7+1.6 
0-1 _2,i 


Fe^t 
Ni«t 


(Zq) 

(Zq) 


0.31 ± 0.02 

1 r,o + 0.50 

^■^'-'-0.91 


0.32 ± 0.03 

-I 08+0-52 



^ from Fe xxvi Ka/Fe xxv Ka. 

^ from (Ni xxvii K«+Fe xxv K/3)/Fe xxvi Ka. 

from Ni xxvii Ka/Fe xxvi Ka. 

from (Fe xxvi Ka-|-xxv Ko)/continuum. 
° from Ni xxvil Ka/continuum. 
t assuming T = 10 keV. 

the higher temperature determined by the broad-band 
spectrum and the much stronger Fe xxv K/3 line. 

As discussed previously, the 2T analysis of Chan- 
dra data suggested that another spectral component is 
needed if no absorption edge modeling is applied. Fig.lS] 
shows the line ratios predicted by the best-flt models 
from J 3.2 over a wide range in temperature for the hot 
phase. Since the spectral energy range used in this flt- 
ting is far enough from the Si edge, it is not necessary 
to modify the spectral model even if the Si edge indeed 
needs to be corrected. The broad-band Chandra spec- 
trum is not sensitive to the hot phase temperature of the 
2T model once it exceeds 15 keV (Fig.|4|. With the good 
constraint from the Fe xxvi Ka/Fe xxv Ka line ratio, 
models with Thot > 20 keV, which are composed of great 
amounts of cooler gas, are rejected. Meanwhile, the ra- 
tio of higher energy states (Ni xxviii Ka, Fe xxvi K/3, 
Fe xxv K7, K5) to the well-measured Fe xxvi Ka line 
suggests that models with lower Thot are preferable. 

As for the 2T models based on Chandra with an ab- 
sorption edge model and XMM-Newton broad-band spec- 
tra, predicted line ratios all agree with the observed 
value. In fact, models with Thot > 20 keV from XMM- 
Newton data are essentially a one temperature model, 
since the normalization of the hot component in these 
models is zero. Adding the fact that an additional tem- 
perature component does not signiflcantly improve the 
of the flt for those spectra and the remarkably good 
agreement on the temperature measured by the contin- 
uum and the iron lines from both Chandra and XMM- 
Newton, we conclude that the simple IT model is ad- 
equate to describe the X-ray emission from the central 
3' region of A1689. 



dances found in S3.1 for Chandra data are likely due to 



4. DEPROJECTION ANALYSIS 

Assuming that the hotter phase gas has the 3D temper- 
ature proflle of L08, the radial distribution of the cooler 
gas can be derived. We extracted spectra from concentric 
annuli up to 8.8' (1.2 /i~^Mpc). The emission from each 
shell in three-dimensional space was modeled with an 
absorbed two-temperature APEC model with Thot fixed 
at the value of |L08| and then projected by the PRO- 
JCT model in XSPEC. Because of the complexity of this 
model, we used coarser annular bins than those used in 
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0.01 




Fig. 6. — The 4.5-9.5 keV Chandra spectrum of the central 
3' region. The spectrum is modeled with an absorbed th erm al 
bremsstrahlung plus the seven Gaussian lines listed in Table[374] 



'L08' Data of'L08'werc binned using the weighting scheme 
of Mazzotta et al.^ (2004) to produce a spectroscopic-hke 
temperature. J cool, abundance, and the normaHzation 
of both components were free to vary. The outermost 
two annuh were background dominated, so spectra were 
binned to have at least 15 net counts per bin at r = 
4.8'-6.5' (625-852 h-^kpc) and 2 net counts at r = 6.5'- 
8.8' (852-1161 h-^kpc) (see |L08| predicted the gas 

temperature only up to 721 /i^^kpc, and that tempera- 
ture was slightly below the observed one. Therefore, we 
allowed Thot to change in the last two bins. The cold 
component was removed and the abundance was fixed at 
0.2 solar in these regions in order to constrain the rest of 
the parameters better. 

Assuming two phases in pressure equilibrium, the vol- 
ume filling fraction of the ith component can be obtained 
from 

_ Norm.i;^ 
•'^ E.Norm.Tj 



(e.g., Sanders & Fabian [2002 1 . Once fi is determined, 
the gas density Pgi = ^e^mpUei can be derived from 



Norm, 



10 



-14 



At:{{1 + z)DaY 



neiiiHifidV, 



(2) 



where nn/ne and /Xg are calculated from a fully ionized 
plasma with the measured abunda nce (He abundance 
is pri mordial, and others are from Anders & Grevesse 
pSgl). For Z = 0.3 Zq, nn/rie = 0.852 and /Xe = 1-146. 
Tig. [9] shows the results of this deprojected 2T analysis. 
The IT modeling, in which emission from each shell has 
only one component, and the results from L08, are also 
shown. If the cluster has a temperature profile of |L08| 
then 70-90% of the space within 250 ft.~^kpc is occupied 
by the "cool" component with a temperature of ^ 10 
keV, based on Chandra absorption edge corrected data, 
a nd this gas constitutes 9 0% of the total gas mass 



Kawahara et al. (20071 show that local density and 
temperature mhomogeneities do not correlate with each 
other in simulated clusters, which undermines the as- 
sumption of two phases in thermal pressure equilibrium. 
However, other cosmological simulations find that gas 



TABLE 8 
Total mass profile 



r 

(/i-ikpc) 


A^lT+cdgG 

(IOI^/i-iMq) 


A^2T+cdgc 
(IOI4/1-IM0) 


M2T 

(IQI-I/i-IMq) 


09+21 

9411? 
2641158 

559 

855+lfg 


"■""-0.01 
n 44+001 

"■^^-0.01 
1 qfi+004 


n nq+009 

u.uy_o 08 

n 67+0-13 

9 fi7+0.60 
^-D'_0.50 
c 90+1.67 
•J X 3 1 

7 qq+4:84 

' -'J'J— 4.27 


0-631°:°^ 
q nQ+0.06 
■^-'-'■^-0.07 
c c-7+0.53 
''■'^'-0.50 

7 nn+3-88 



Note. — 2T assumption is only held within 625 h^^kpc. 
The upper and lower limits of r indicate the radii f of two 
contiguous rings used to calculate the mass. See text for defi- 
nitions of r and f. 



motions contribute about 5-20% of t 


he total pressure 


support (e.g.. 


Faltenbacher et al. 12005 


Rasia et al. 12006 


Lau et al.|l2009'). It the pressure balance is ott by 20%, it 



or the volume filling fraction (< 8%). 

5. MASS PROFILE 

Given the 3D gas density and temperature profiles, the 
total cluster mass within a radius r can be estimated 



from the hydrostatic equilibrium equation (e.g., Sarazin 
19881, 



M{r) 



kT(r)r fdlnpgir) dlnTir) 
Gji nip \ d In r d In r 



(3) 



For Z = 0.3 Zq, m = 0.596. 

If the gas has two temperatures with two phases in 
pressure equilibrium, the total mass still can be derived 
from Eq. [3] with pg , T replaced by Pg^^^^ , Thot , respec- 
tively. 

5.1. Nonparametric method 

To evaluate the derivatives in Eq. |3] we took the dif- 
ferences of deprojected temperature and the gas density 
in log space. The radius of each annulus was assigned at 
f such that 



F3D{r) — 



'3 _ ^3 ^ 

V out in J 



F3D{r)i7Tr^dr, (4) 



where F^jj is the deprojected flux density from a finely 
binned surface brightness profile, and r^^ (rout) is the 
inner (outer) radius of the annulus. The radius r outside 
of the brackets of Eq. [3] is taken at the geometric mean 
(i.e. the arithmetic mean in log scale) of the radii of 
two adjacent rings, r = ^fi f^+i , and the temperature 
is linearly interpolated at this radius. Because errors 
from e.g. T and dT jdr are not independent, standard 
error propagation is not easily applied. Uncertainties 
are estimated from the distribution of 1000 Monte-Carlo 
simulations of T and pg profiles. Fig. [9] shows the total 
mass profile from both IT and 2T models and the results 
arc listed in Table|8] Two-temperature modeling, based 
on the Thot of L08, increases the total mass by 30-50% for 
all radii within 625 /i~^kpc. Beyond that radius, the 2T 
assumption is not held because of the lack of constraint 
on Thot- 

Although the inclusion of an absorption edge in the 
spectral model greatly changes the derived composition 
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Fig. 7. — The predicted IT plasma line ratio (dotted line) as a function of temperature, for various lines. The observed ratio and its Icr 
confidence are shown as a solid line and shaded region. The circles show the fitted results of 100 simulated Chandra spectra drawn from a 
VAPEC model with 9 Ni© /Fe© . 
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The solid line and shaded region shows the observed ratio and its la error, inc x-axis is in log scale. 



of the multi-phase plasma, it does not affect the mass 
measurement much. This is because we use a fixed Thot 
profile. Once the temperature is determined, the total 
mass only depends on the logarithmic scale of the gas 
density, which produces ~ 13% difference at most. 

5.2. IT parametric method 

If the temperature does not vary dramatically on small 
scales, we can obtain a mass profile with higher spatial 
resolution since the gas density can be measured in de- 
tail from the X-ray surface brightness with the assump- 
tion of a certain geometry of the cluster. To achieve 
this, modeling of the temperature and the gas density 
is nece ssary. Following the procedure of Vikhlinin et al. 
( 2006 1 , we project the 3D temperature and the gas den- 



sity models along the line of sight and fit with the ob- 
served projected temperature and th e surface b rightness 
profiles. A weigh ting method by ,Mazzotta et al. (2004j , 



Vikhlinin ( 2006 1 is used to predict a single-temperature 
tit to the projected multi-temperatur e emission from 3D 
space. This method has been shown ( Nagai et al.||2007 l 

to accurately reproduce density and temperature profiles 
of simulated clusters. 
The gas density model is given by 



{r/rc 



(l + r2/r2)3'3-"/2 (l + r^/r^-yy/'y 
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Fig. 9. — Temperature, emission measure ratio of the cool component EM^^^i/ EMtot, volume filling fraction of the cool component 
/cooli gas number density rie, cumulative gas mass Afgas, and cumulative total mass M^o profiles from the 2T deprojection analysis with 
an absorption edge correction (squares) and without the correction (circles). Also shown is the IT analysis (diamonds) and results from 
Lemze et al. (2008a) (asterisks). Tj^ot of the first 4 annuli was fixed at the value derived from lensing and X-ray brightness data ( Lemze] 
elT al.^2UU8a ), which were grouped into fewer bins. The cool component of the last 2 bins was frozen at zero. The 2T assumption is hela 
within 625 /i~^kpc. X-data points of the IT and 2T models have been shifted by +10% and -10% for clarity, and their error bars are also 
omitted. 



which orig inates from a (3 model fCavaliere & Fusco- 
Femiano 1978) modi fied by a power-law cus p and a steep- 
ening at large radii (Vikhlinin et al. 1999). The second 
term describes a possible component in the center, espe- 
cially for clusters with small core radius. The tempera- 
ture model is given by 



Tauir) = To 



Dol) 

+ (r/r. 



cool J 



which is a broken power law with central cooling (Allen 
et al.|[200T ). Best-fit parameters for thegas density and 
temperature profiles are listed in Tables [9]and [Tol respec- 
tively. Errors are estimated from the distribution of the 
fitted parameters of 1000 simulated projected tempera- 
ture and surface brightness profiles generated according 
to the observed data and their measurement uncertain- 
ties. Since parameters are highly degenerate, some of the 
best-fit values are not covered by the upper or lower lim- 
its with the quoted confidence level (upper/lower bounds 
are for one parameter). The observed temperature and 
surface brightness profiles, the best-fit model, and the 
surface brightness residual are shown in Fig. |10[ The 
model describes the data very well (x^ /do f =154.3/155). 
The best-fit T3D and Ue models are shown in Fig.jTT] Also 



plotted are the profiles from the spectral deprojection 
fitting (^5.1). Compared to this nonparametric result, 
modeling 1 3D and can avoid fiucutations from the di- 
rect spectral deprojection, which is a common problem 
as the deprojecti on tends to amplify the no ise in the data 
(see Appendix in Sanders fc Fabian||2007 1 . 

Although the second break ot the first term in Eq. |5] 
was desig ned to describe the steepening at V g > 0.3r2oo 
( fVikhlinin et al.l|1999[ |2006[ |Neumann|(2005[ ), we found 
that it the initial guess tor is not big enough, Ts tends 
to converge to a relatively small value, » 200 h~^kpc, 
compared to the typic al value of 400-3000 h ~^kpc for 
nearby relaxed clusters (IVikhli nin et al.|2006 l. It is pos- 
sible to use the first core radius Tc or the core radius of 
the second component to account for the sharpening 
at 200 h~^kpc. This consequently yields a more reason- 
able Ts at « 1 h^^Mpc. Both cases, small (Model 1) 
and large (Model 2) r^, give acceptable fits with x^/dof 
of 153.4/155 and 154.3/155, respectively. However, large 
is harder to constrain. This makes the mass estimate 
more uncertain at large radii than the small case. 
Comparing the surface brightness profile of t he north- 



eastern (NE) part t o the southwestern (SW), Riemer 



S0rensen et ah] p009) found that the NE part is 5-15% 



brighter outside 350 h ^kpc and 25% under-luminous at 
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Fig. 10. — Projected temperature and surface brightness profiles 
with the best-fit model (solid lines) and its 90% confidence bounds 
(dashed lines). Bottom panel: residual between the surface bright- 
ness and the model. This fit gives a jdof of 154.3/155. 
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Fig. 11. — Best-fit T3D and ne models (solid fines) and 90% 
confidence bounds (dashe d lin es) . Also shown arc unparameterized 
results (diamonds) from i]5.1[ 



70 /i~^kpc than the SW. To see if this asymmetry can 
affect the mass estimate, we fit a symmetric model to 
the image and iteratively removed any part of the clus- 
ter that deviates significantly from the azimuthal mean, 
mainly the northern clump at 460 /i~^kpc, the south- 
ern less luminous region at 330 /i~^kpc, and possibly 
some point sources not completely removed beforehand. 
We did not exclude these regions from our temperature 
measurement since they were unlikely to bias the aver- 
age temperature much for such a hot cluster, as shown in 
Fig. [4] that at least 10-20% of the total emission measure 
from another spectral component was needed in order the 
change the spectroscopic temperature by 1 keV. Best-fit 
gas density and temperature for these models are listed 
in Tables [9] and 10 labelled with Model 3 (small r^) and 
4 (large rT). 

The total mass profiles from these analytic gas density 
and temperature models are given in Tablefll] We list 
the total mass at the radii where masses from the non- 
parametric method are evaluated (Tablelsl. The last en- 
try of TablefTT] shows the total mass at the boundary of 
the ACIS-I cliips, 12' (1.6 /i^^Mpc « raoo), where Sx is 
detected at < Icr. Removing asymmetric parts from the 
image or restricting to be greater than 350 /i~^kpc in- 
creases the total mass estimate with < 10%. Nonetheless, 
these differences are not significant. We combine samples 
of the best-fit parameters of substructures removed cases 
(Model 3 and 4) as our best result. 

Fig. [12] shows a comparison of parametric and non- 
parametric mass profiles. The non-parametric mass pro- 
file is from a more finely binned deprojected data than 
those shown in Fig. [9] Results from these two methods 
are fully consistent with each other, although their er- 
rors are quite different. For the non-parametric method, 
we simply assign the observable, e.g. dT/dr, at a cer- 
tain radius, so the uncertainty associated with the posi- 
tion is not included in the error on the mass, ctm, but 
separately shown on the radius. Therefore, ctm appears 
smaller if data are binned more coarsely. For the para- 
metric method, the dependency of ctm on the data bin- 
ning is weaker. The departure from the model for any 
data point is assumed to be random noise and is filtered 
out through the fitting. Hence, au reflects only the un- 
certainty of the fitted function and it depends strongly 
on the modeling. 

5.3. Comparison with other studies 
T he to tal mass profiles of |A04[ based on XMM-Newtc 



and LOS a joint X- ray , strong and weak lensing study are 



also show n in Fig. 12 Our r esult is in good agreement 
with A04 but disagrees with LOS around ~ 200 h~^kpc. 
To compare our mass estimate with other lensing works, 
we derived the total mass density and integrated it along 
the line-of-sight. The total mass density, p, is obtained 
through the hydrostatic equation, 

k 



Pg-rrip 



{\I'T + TV^lnpg + Vlnpg • VT) . (7) 



For the nonparametric method, we evaluated Eq. [7] in 
a similar fashion as we did in j ]5.1[ Errors were esti- 
mated from the Monte-Carlo siniulations of deprojected 
T and pg profiles. Fig. 13 shows the surface mass density 
profiles from both parametric and nonparametric meth- 
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TABLE 9 

Best-fit parameters for the gas density (Eq[5]i 



no Tc Ts 


Q 


/3 


e 


no2 r-c2 


/32 


7 


IQ-'^hhcra-^ lO^ft-ikpc lO^fe-ikpc 








10"l/i3cm-3 lO^/i-^kpc 







n~l -^fiS+O-Ol nQ^+O-lO 1 QQ+O-19 r, 07+0.32 n -^fi+O-O" 1SQ+0-24 n Qn+O'^S 9 n4+''-26 yin+^-S* 4 fi/1+1'23 

(i) cl.b8_o 7g U.^&_o oo i.yy_o U.dl'.g.ig U.db_Q 03 i-Sy_o.28 U.^U.g.oa i.U4_(, 49 '■^^>-2.63 4.b4_(, gg 

C9~l n 1 fit;+0-ll 11 1 1/l+0-08 70+0.02 n c;n+3-66 n 97+0.03 ,79+0.16 fiin+0-85 9 t;fi+7.06 

(^j U.f-S.Qoy ^-03-0.21 iJ^-O+0.1 ^■^'^-O.IQ U.('d_oQ4 U.&U_o4g i./'^_g72 0-^^-3.85 ^■30-0.28 

("?~) 2 46+1-°** n4fi+°-°3 2 7n+°-l^ (,70+0.09 „ 44+0.02 1 Q1+0.36 (,,q+0.02 „ 10 0+°-°'^ 7 52+0*6 

(A) ^-^0-0.22 "-^0-0.15 ^-'^-0.15 '-'■'''-0.39 "'^^-O.OS ^^-"^^-0.09 '-'■^^^-0.04 ^•'^'^-0.89 ^^"•'-'-6.1 '■'^^-2.56 

IA\ nfi4+0-33 1 QS+0-80 ip, 9+27.9 +0.11 +0.19 , .^+2.91 ^ 97+O.O9 ,77-0.19 9^0+8.62 

(4) U.b4_g i.y8_o.6o J^U.^_6.72 ^-^^-O.bd ^-"-0.12 ^-40-1.33 ^-^'-0.02 ^■"-1.22 O-Ui_5.08 ^-^8-0.71 

Note. — (1) small rg, (2) large Vs. (3) small Ts with northern clumps removed, (4) large Ts with northern clumps removed. Errors are 95% CL 
for one parameter from 1000 Monte-Carlo simulations. Since parameters are highly degenerate, some of the best-fit values arc not covered by the 
upper and lower limits at this confidence level. 
^ parameters hit the hard limit. 



TABLE 10 

Best-fit parameters for the temperature (Eq[6| 



To 


Tmiii/To Tcool n 


a 


b 


c 


d 


kcV 


lO^/i-ikpc lO^ft-^^kpc 











(1) 12.7! 
(2) 

(3) 14.4 

(4) 11.8 



+6.9 
2.8 

12 1+'^'^ 
+5.8 
5.3 
+8.4 
2.9 



,+0.25 
-0.34 
7+0.11 
-0.41 
hO.60 
-0.13 

0.86l°o:ro 



0.73: 
0.87j 
0.40^ 



0.67 
0.64 
0.77 
0.74 



1.66 

0.29 
1.75 
0.15 
1.42 
0.28 
+ 1.84 
0.29 



12.5 
8.36 
28.6 
9.08 



+20.8 

-6.8 
+21.7 
0.93 
+0.5 

25.0 
+23.1 
2.42 



0.02 
-0.02 

0.11 
-0.01 



+0.11 

0.14 
+0.16 
0.11 
+0.10 
0.24 
+0.17 
0.12 



0.86 
1.37 
0.40 
1.68 



+0.64 
0.36 
+0.65 
0.80 
+1.11 
0.03 
+0.81 
1.15 



0.81 
0.89 
0.59 
1.32 



1.82 
0.40 

+2.11 
0.43 
1.56 
0.19 

+ 1.64 
0.92 



2.96 
6.78 
1.44 
7.07 



+2.14 
2.45 
0.83 
6.30 
+3.86 
1.17 
0.61 
6.57 



Note. — (1) small r^, (2) large r^, (3) small with northern clumps removed, (4) large with northern 
clumps removed. Errors arc 95% CL for one parameter from 1000 Monte-Carlo simulations. 



TABLE 11 
Parametric total mass profile 



r {h -"^kpc) 




M{r) (10"/i- 


-IM0) 






(1) 


(2) 


(3) 


(4) 



32 

94 
264 
559 
855 
1579 



0.07 
0.46 
2.01 
4.15 
5.73 
8.60 



+0.01 
0.01 
+0.01 
0.02 

+0.07 
0.03 
+0.32 
0.02 
+0.46 
0.05 
+0.59 
1.02 



0.07 

0.46 
1.92 
4.43 
6.37 
9.52 



0.01 
0.01 
+0.01 

0.03 
+0.09 
0.02 
+0.18 
0.11 
+0.35 
0.36 
+2.22 
-0.97 



0.07 
0.47 
1.92 
4.64 
6.41 
9.85 



+0.01 
0.01 
+0.02 
0.02 

+0.09 
0.03 
+0.27 
0.03 
+0.40 
0.23 
+0.56 
1.52 



0.07 
0.44 
1 

4.86 
7.03 
9.55 



+0.01 
0.01 

+0.03 
0.02 
+0.08 
-0.08 
+0.15 
0.51 
+0.71 
1.00 
+3.31 
1.43 



Note. — (1) small r^, (2) large Vs, (3) small with north- 
ern clumps removed, (4) large rs with northern clumps re- 
moved. Errors are 68% CL from 1000 Monte-Carlo simula- 
tions. 

od s, along with the HST/A CS strong lensing analysis 
of iBroadhurst et al. (2005a I, an d the combined Subaru 
distortion and depletion data by [Umetsu fc Broadhurst| 
( |2008J . Since it requires at least 3 points to calculate the 
second derivative, p at the boundary is unknown. This 
will introduce additional systematic errors to the inner 
and the outer projected profile. To demonstrate how this 
may affect our nonparametric result, we insert two arti- 
ficial points at 1" (2 /i^^kpc) and 13' (1.7 h~^Mpc) to 
the nonparametric T and pg profiles with their values es- 
timated from the parametric model. The projected den- 
sity derived this way is shown in red filled diamonds in 
Fig. 13 The X-ray data are consistent with those from 
the weak lensing, but disagree with the strong lensing 
analysis. Although the nonparametric data appears to 
agree with the strong lensing estimate at r = 80 h~^kpc, 
this is probably due to the temperature fluctuation men- 
tioned in |5.2[ 
The mass discrepancy is manifested when comparing 
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Fig. 12. — The parametric mass profile (solid line) compared to 
the unparameterized result (diamonds). Dashed lines show the 95% 
confidence bounds. Also shown are XMM-Newton result from lA0'4l 
(cros ses) and combined X-ray, stron g and weak lensing analysis of 
|L08| (asterisks). The mass profile of |L08| is mainly determined by 
the lensing data. 



the cumulative projected mass profiles, M2 D, shown in 
Fig.[T 4| Th e weak lensing M2D profile of |Umetsu~Sg] 
Bro adnurst ([200 8) includes the integration of the data 
of |t}roadhurst et al.| ( |2005a| in the inner region. Un- 
certainties are from Monte-Carlo simulations of the con- 
vergence pr o files. The last 3 data points of [Umetsu IE] 
Broadhurst (20081 (1-2.3 /i~^Mpc) are discarded since 
only the upper limits are available. Also shown are 



parametric strong len sing profiles (Halkola et al. 2006 



Limousin et al. 20071, and other X-ray analyses (A04; 
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Fig. 13. — Surface mass density profiles from non-parametric 
(open diamonds) and parametric X-ray model (solid and dashed 
lines, 95% CL), compare d to HST/ACS strong lensing analysis of 
[Broadhurst et ar| | |2005a| (triangles), and combined Subaru distor- 
tion and depletion data by|Umetsu & Broadhurst (2008), based on 
a maximum entropy method (circles). t"'illcd diamonds show the 
mass from the nonparametric T-jo and rie profiles that include es- 
timations from the parametric result at 1" (2 /i~^kpc) and 13' (1.7 
/i-^Mpc). 



Riemer-S0rensen et al.||2009l. To convert M^d to M2D, 
assume that the lastTaata point reached the clus- 
ter mass limit, which unavoidabl y leads to underestima- 
tions especially at large radii. |Riemer-S0rensen et al. 
( 2009 1 use only the SE part of the cluster and four of the 
Chandra observations (excluding ObsID 540) and derive 
M2D based on a best-fit NFW model fit to the M^d 
profile. Their mass profile is generally lower than our 
estimate at most radii. This is contradictory to most 
findings that claim that the hydros tatic mass is underes 



timated in unrelaxed systems (e.g., Jeltema et alT][2008 1 . 
Using such reasoning, and removing the INli; part, pre- 
suma bly disturbed according to Riemer-S0rensen et al.| 
(|2009| , should increase the overall mass estimate. The 
A-ray M2D is 25-40% lower than that of lensing within 
200 /i-^kpc, corresponding to a ~ 1.4 x 10" h'^MQ dif- 
ference in the total projected mass. 

5.4. NFW profile parameters 
The total mass prof ile M^jj was fit to the NFW model 



(Navarro et al. 19971 to obtain the mass and the con- 
centration parameter. To fit the nonparametric data, we 
weighted each point according to its vertical and hori- 
zontal errors, given by. 



2 fdMV 



(8) 



where ar is assigned to be 68% of the width of the hori- 
zontal error bar and dM/dr is iteratively evaluated from 
the NFW model until it converges. In the parametric ap- 
proach, a NFW model was fit to the parametrized mass 
profile that evaluated only at the radii where the pro- 
jected temperature was measured with errors estimated 
from the standard deviation of a sample of mass profiles 
constructed from the simulated Tpioj and Sx profiles de- 



Q 
CM 



0.1 

I I 1 1 
- X 

□ 
▲ 



r (arcmin) 
1 



This wo* 

Andersson & Madejsl(i 04 
Riemer-Sorensen et al. 08 
Broadhurst et al 05 
Umetsu & Broadhurst 08 
Halkola et al. 06 
Limousin et al. 07 




100 
r (/7-'kpc) 

Fig. 14. — Projected mass profiles from non-parametric (open 
diamonds) and parametric analyses (solid and dashed lines, 95% 
CL), compar ed to XMM-Newton result from A04 (crosses), Chan- 
dra result by 'Riemer-Sorensen et al. (2009) (squares), HST/ACS 
and Sub aru results by Broadhurst et al. (2005a) (triangles) and 
|Umetsu~fc Broadhurst (2008) (circl es). We integrated the lensing 
surtace mass profile (shown in Fig. |13| ) and estimated its uncer- 
tainties from Monte-Carlo simulations. Also shown are parametric 
strong lensing profiles of Halkola et al. ( 2006 ) and Limousin et aT] 
l |2007| l (shaded regions, 68% CL). Kiemer-S0rensen et al.^ (2003| l 
used only SW part of the X-ray data and converted M^o to M2r> 
with a NFW profile. |A04| assumed that the last data poin t reached 
the cluster mass limitTTilled diamonds, same as Fig.|13[ 



mass profiles in the sample. Resulting NFW parameters 
were used to estimate the uncertainty. 

Table 12 lists the best-fit NFW parameters, M200 and 
C20O1 for the total mass from both methods and from 
other studies, all converted to the adopted cosmology. 
Compared to other X-ray studies, our derived M200 is 
30-50% higher, closer to weak lensing results. The differ- 
ences between our NFW parameters and those of |A04| 
from XMM-Newton are primarily attributed to their 
slightly lower but yet consistent mass at the last data 
point (Fig. 12). This demonstrates that the accurate 
mass measurement at large radii, where systematic errors 
are usually the greatest, is crucial to the determination 
of NFW parameters. 

Our results are consistent with weak lensing measure- 
ments, but with a lower concentration than what recent 



hurst 


2008 


are ac 


Ided \ 



JCorless et al. 2009). When these analyses 
are added with strong lensmg information, a very tight 
constraint on the concentration par ameter can be ob- 
tained, giving 

C200 = 9.9lE]j (Umetsu & Broadhurstl 
which hardly can be reconciled with our value, 
5.3^J'2. However, if the gas emission i s mo deled with 
two spectral components with Thot from |L08l the X-ray 
derived concentration is in a closer agreement to those of 
combined strong and weak lensing studies, but this also 
implies that the majority of the gas is in the cool phase 
and occupies most of the intracluster space (j |3.2| . 

5.5. Gas mass fraction 
The cumulative gas fraction /gas — -/Wgas/Mtotai, de- 



scribed in ^5.2 We repeated this procedure for all of the rived from our best-fit T3D and rie model, is 0.098 



+0.003 
0.004 
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TABLE 12 

Comparison of best-Fit NFW Parameters 



Method 


Instrument 


A/200 


C200 




Reference 






(10^5^-1^0) 









Spherical model 



X-ray (IT+edge) 


Chandra 


^ ^ ^-(-0 45 




6.3/8 


X^-ray (parametrized Tsq, ?2e) 


Chandra 


94+0-" 


6 6+°1 

0.4 




X-ray {2T'') 


Chandra 


— U.ZO 


7.6+li 


2.2/3 


Tv-i ay 1 z ± n^eu^e i 


v_>iianQra 


, 19-1-0.53 

^•^^-0.29 


q cj+0.7 


u. 1/ 


X-ray (IT) 


XMM-Newton 


0.63 ± 0.36 


7 6+l'^ 
'•0-2.6 


7.6/8 


X-ray (IT) 


Chandra 


0.55 


10.1 <^ 


1.6/13 


SL 


ACS 


2.29 


6 3+1** 
"•■^-1.6 




SL 


ACS 


2.16 ±0.32 


5.8 ±0.5 


0.8/11 


WL 


CFHT 


0.97 ±0.13 


7.4 ± 1.6 


WL 


CFHT 


0.90 ±0.17 


13.1 ± 7.5 




WL 


Subaru 


1.24 ±0.14 


10.5t^;* 


332/834 


SL-I-WL 


ACS+Subaru 


1.22 ± 0.13 




13.3/20 


SL-I-WL 


ACS+Subaru 


1.31 ± 0.11 




335/846 


SL-l-WL+X-ray {Sx) 


ACS+Subaru+Chandra 


1.42 


g 7+0.8 


15.3/24 


Triaxial model 










SL+WL b 


ACS-l-Subaru 


1 1 C-+0.26 
-L--LO_0.45 


1 4+1-8 


378/362 


WL'^ 


CFHT 


0.83 ±0.16 


12.0 ±6.6 





this work 
this work 
this work 
this work 



A04 



2009 1 



Hiemcr-S0rcnsen et al 
Broadhurst ct al. ( 20 05a| ^ 
Halkola ct al. (i!U06) 
Limousin ct al. (ZUU7| 
Corlcss ct al. ( 20uy ) 
Umctsu & Broadhurstjj 2008 1 
Broadhurst ct al. ( 2005 
Umctsu & Broadhurst 



LOS 



(20081 



|Oguri et al. 


(2005 




Corless et ai 


.|I|2ULI9|| 



Note. — see |Comerford & Natarajan||2'007[ ; |Umetsu & Broadhurst| ( |2b08[ l; [Corless et al.| ( |2009[ l for a more complete compilation. 

^ with Thot from |L08| 

^ under a flat prior on the axis ratios. 

under a prior on the halo orientation that favors the line-of-sight direction. 

converted from best-fit parameters, po = 7.79 X lO^M0kpc~^, = 174 kpc {h = 0.7, fim = 0.28, and = 0.72), of |Riemer-S0rensen| 
|et al.|f2009[ , not consistent with their quoted value of 5.6 since they did not adopted the commonly defined C20O5 the concentration at r2oo 

where the enclosed mean density is 200 times the critical density (private communication). 



-1.5 



ik-jQ at r25oo (493_;^Q /i^ kpc), ^20% higher than what 
is found using XMM-Newton data ( ,A04^ . In spite of this 
seemingly large difference, the data agree that /gas does 
not converge at r25oo- Much Hke in the case of A1689, 
the low-z relaxed cluster A1413 does not have a strong 
cooling core and also has a steadily rising /gas profile 
out to rsoo ( Pratt fc Arnaudn2002j . Comparing the /gas 
profile of A1413 wit h anotEer nearby prominent cooling 



TABLE 13 

Comparison of Afsoo 



core cluster, A478, Pointecouteau et al 
late that the flatter / 



( |2004[ ) specu- 
protile of A478 is related to the 



Method 


Af500 

(lO^/i-iMo) 


'^'500 

(h-iMpc) 


parametrized Tsd, 
M500 - Tx t 
M500 - t " 


7 0+1.3 

7.7 ±0.2 

7 7+0.5 
'•'-1.2 


i-'JJ^-o.os 
1.03 ±0.01 

1 03+°'°^ 



presence of a cooling co re. Our /25qo is 11 % lower than 
the mean gas fraction of Allen et al. ( 2008 ) derived from 
42 relaxed clusters observed with Chandra, but our /500, 
0.12 ± 0.01 /i7n^ ■^ agrees with in 1% of the M - /gas rela- 



Vikhlinin et al. 



120091 with 



tion of Vikhlinin et al. (2009 1 



DISCUSSION 



Nagai ct al. (20071 s how t hat following the data analy- 
sis of |Vikhlinin et aL (20061, the hydrostatic mass is un- 
derestimated by 14 ± d% within estimated r25oo for sim- 
ulated clusters visually classified as " relaxed" . Based on 
the X-ray morphology, A 1689 is likely to be categorized 
as a relaxed cluster. The X-ray centroid is within 3" of 



t Scaling relations from 

indices fixed to self-similar theory values. Errors only 
reflect the measurement uncertainties. Dispersions of 
the relation is not included. Tx = 10.1 ± 0.2 keV, 
measured from r = 1. 14' — 7.6' (~ O.lSrsno — ''500 )• 

By solving Eq. 14 of | Vikhlinin et al. | pOOg] ! . The fi- 
nal Yx = Tx X Afgas determined at rsoo is (5.ll5'2) ^ 
1(^14 /i-2.5]Vlg ]^-gV. 



unusual for "unrelaxed" clusters, referring to those with 
secondary maxima, filamentary structures, or significant 
isophotal centroid shifts. 



Table 



13] shows the comparison of measured M500 with 



others derived from the M. 



the le nsing and optical centers (Andersson & Madejski 
2004 1 , with a very min imal centroid shift or asymmetry 



Yx and M, 



500 



zUU4 1 , With a very min imal centroid shiit or asymmetry 
( [Hashimoto et al.|2007 1 . At the X-ray estimated r25oo of 
493 /i~^kpc, we derive an enclosed hydrostatic mass of 
(4.2 ± 0.3) xK)"/i-iM0, « 30% lower than the lensing 
mass from|L08!. At r = 200 /i~^kpc, this becomes a 50% 



Tx rela- 
tions of Vikhlinin et al. ( 2009 1 , calibrated from 49 low- 
z and 3y high-z with (z) = 0.5 clusters observed with 
Chandra and ROSAT. A very good agreement has been 
achieved between these estimates. Since the M500 — Yx 
relation is insensitive to whether the cluster is relaxed 



difference (see Fig. 12 1. Such a strong bia s is no t seen in 
the relaxed cluster sample of iNagai et al. (20071, assum- 
ing that the lensing mass is unbiased, although this is not 



or not (Kravtsov et al. 2006 ) and merging clusters tend 
to be cool for their mass ( jMathiesen fc Evrard 20011, 
consistency among these rriass estimates indicates that 
A1689 is relaxed in the sense that it behaves like other 
"relaxed" clusters on the scaling relation. 
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On the other hand, projection effects, such as triaxial 
halos or chance ahgnments, always have to be taken into 
account when comparing projected (lensing) and three- 
dimensional (X-ray) mass esti mates. F rom kinematics of 
about 200 galaxies in A1689, Lokas et al.| (j2006) sug- 



gest that there could be a few distant, possibly non 
intera cting, substru c tures s uperposed along the line of 



sight. Lemze et al. (2008b I, based on a 0.5 x . 5 deg 



VLT/V1M(JS spectroscopic survey from |Czoske| ( |2004 1 

which includes ^ 500 cluster members, disagree with this 
projection view. They conclude that only one identifi- 
able substructure at +3000 km/s, 1.5' to the NE (The 
X-ray clump is at 3.5' NE). This backgr ound group 



is seen in th e strong lensing mass analysis (Broadhurst 
et al. 2005a I, but is determined not to be massive (< 10% 



of the total mass in the strong lensing region). Nonethe- 
less, the higher than usual velocity dispersion in the clus- 
ter center, ~ 2100 km/s, indica tes that the central part 



is quite complex (Czoske 



2004). This may also imply 



that the halo is elongated in the line-of-sight direction 
as galaxies move faster along the major axis. 

For powerful strong lens systems, like A1689 
sphericity i s never a justified assumptio n (e.g., 
etaL][2007| . |Oguri fc Blandford| ( [2009] ) show t 



halo 



Hennawi 
tiat these 

superlens'' clusters almost always have their major axes 
aligned along the line-of-sight, with more circular ap- 
pearances in projection and ~ 40 — 60% larger concen- 
tratio ns than other cl usters with similar masses and red- 
shifts. Gavazzi (2005) demonstrates that using a prolate 
halo with axis ratio 0.4, they were able to explain 
the mass discrepancy between the lensing and X-ray es- 
timates of cluster MS 2137-23. This cluster has a well 
defined cool core (e.g., Andersson et al.|[2009 l, thus pre- 
sumably relaxed, and yet a factor of 2 ditterence in the 
mass is not lessened with a multiphase model for the core 
region (Arabadjis et al. 2004). In contrast, triaxial mod- 
eling not only solves the mass inconsistency, but also the 
high concentration problem and the misalignment be- 
tween stellar and dark matter components in MS2137-23 
( [Gavazzi 2005). 

'ib see how the triaxiality changes our mass measure- 
ments, we modeled T3D and Pg with prolate profiles, by 

replacing r in Eq.jHjandjejwith (x^/a^+y^/fe^-l-z^/c^)^/^, 
where we assumea a = b < c and the major axis, z- 
axis, is perfectly aligned along the l ine- of-sight. Follow- 
ing the same analysis outlined in S]5.2| but with differ- 



ent projection factors, we obtained best-fit Tsd and pg 
profiles. The derived mass prof iles under various axis 
ratios a/c are shown in Fig. 15 The uncertainties on 
S(r) and M2D{'r) are similar to those in Figs. 
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13 and 



We integrated the density from z = — 4.5/i^ Mpc to 
5ft-^^Mpc (« 3r2oo for a/c=l) for all the cases. The 
uncertainties of T3D and pg profiles at large radii (> 10% 
aX r — rsoo and increasing further afterward) does not 
significantly change the projected mass at smaller radii 
(< 3% within 500/i^ikpc). 

The total mass enclosed within a sphere of radius r, 
M^£){r), and the spherically averaged mass density, p{r), 
are basically unchanged under different assumptions of 
triaxiality, considering the typical measureme nt uncer- 
The sam e conclu s ion w as drawn by IPiffaretti 



tainty . 



et al. 
7?~br a 



(20031 and Gavazzi (2005), though they assumed 



WFW model with gas isothermality. For the 



r (arcmin) 
0.1 1 




r (h^'kpc) 

Fig. 15. — Best-fit mass profiles for various axis ratios a/c from 
Model 1 (dash-dot line), 2 (dotted line), 3 (solid line), and 4 
(dashed line) . Top left: total mass enclose within a sphere of radius 
r, A/3j3(r). Top right: spherically average mass density p{r). Bot- 
tom left: azimuthally average surface mass density S(r'). Bottom 
right: projected mass within a cylinder of radius r, -Af2i3(r). 



azimuthally averaged surface mass density ^^(r) or the 
projected mass within a cylinder of radius r, M2D{f), 
a factor of 2 or more difference can be easily made by 
increasing the ellipticity. An axis ratio of 0.6, giving 
M2d{< 45") = 1.4 X lO"/i-iM0 (by a factor of 1.6 in- 
crease), can resolve the central mass discrepancy, but 
overpredicts the mass by ~ 40% at large radii. For a 
ratio of 0.7, the X-ray mass estimate data agrees with 
those of strong and weak lensing within 1% (— Icr) and 
25% (-t-lfj), respectively. Since the gas distribution is 
rounder than that of t he DM , a larger axis ratio than 
the finding of Gavazzi (2005) is expected. 

Not only does the projected mass increases with the 
triaxiality, but also does the steepness of the profile. 
This explains a higher than X-r ay derived concentra- 
tion from the lensing data (^5.4). Although some at- 
tempts have been made to model the lensing mass profile 
with a 3D triaxial halo (Oguri "erar[[2005| ' 



Corless et al. 



2009 ) , no significant constraint on the concentration pa- 
rameter is obtained (Table 12). To break the degeneracy 
between the triaxiality and the concentration, observa- 
tions from different prospective projections, such as X- 
ray, Sunyaev-Zel'dovich effect, or galaxy kinematics, are 
al ways needed. 



Riemer-S0rensen et al. 



^ ([2009f conclude that A1689 

harbors a cool core based on the radial temperature pro- 
file and a hardness ratio map. From this, they conclude 
that A1689 is relaxed, excluding the NE half of the clus- 
ter, where there is a low mass substructure. Based on the 
derived temperature profile we disagree that A1689 con- 
tains a cool core. The temperature of the cluster varies 
radially from 9 to 11 keV with a slight drop only af- 
ter 500 /i~^kpc. This can not be characterized as the 
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properties of a cool co re cluster. In fact, as shown in 
Andersson et al. (20091, A1689 is an intermediate stage 
cluster in terms of central baryon concentration with a 
minimal core temperature drop. This, does not neces- 
sarily provide evidence that the cluster is disturbed but 
we do not cither expect the properties of a cool core 
cluster. Hardness-ratio maps are very sensitive to ac- 
curate background subtraction, especially for high en- 
ergy splittings. We suspect that the ha rdness ratio map 
{S/H = ^[0 3 - Q.Q\/E[Q.Q - 10.0]) in [Riemer-S0rensen 
et al. (2009), used as evidence for a cool core, sutters 
Irom residual background. The ratio decreases rapidly 
with radius from a central value of 2.2 and reaches 0.4 
already at 3'. This is an extremely low ratio for any 
reasonable cluster temperature and it is in disagreement 
with the observed temperature profile. For comparison, 
a background-free spectrum from an isothermal cluster 
at 10 keV would exhibit a count-ratio of ~ 47 in ACIS-I 
given the energy bands mentioned above. The usage of 
unsubtracted hardness-ratios in these bands shows that 
the high-energy band has a significant fractional back- 
ground contribution and hence, is more spatially flat 
compared to the low energy band. This does not provide 
information about the spatial distribution of gas temper- 
atures in the ICM. 

7. SUMMARY 

We have investigated a deep exposure of Abell 1689 
using the ACIS-I instrument aboard the Chandra X-ray 
telescope. In order to study the discrepancy of the grav- 
itational mass from estimates from gravitational lensing, 
to that derived using X-ray data, we test the hypoth- 
esis of multiple temperature components in projection. 
The result of a two-temperature model fit shows that it 
is very important to take into account all details of the 
calibration of the instrument. We detect an additional 
absorption feature at 1.75 keV consistent with an ab- 
sorption edge with an optical depth of 0.13. In analyzing 
multiple additional datasets, we find similar parameter 
values for this edge. 

If the edge is not modeled, fitting the cluster data 
within 3' strongly favors an additional plasma compo- 
nent at a different temperature. However, when this ab- 
sorption feature is modeled, the second component does 
not improve the statistic significantly and the fit results 
is in better agreement with the XMM MOS data. In all 
cases, a second component has to have T > 5 keV in 
order for the hot component to agree with the cluster 



(|2008a|) which is 



temperature predicted by Lemze et al. 
derived from lensing and Sx profiles, 'i'his contradicts 
the assertion that cool clumps are biasing the X-ray tem- 
perature measurements since these substructures would 
not be cool at all. We also find that, if the temperature 
profile of the ambient cluster gas is in fact that of jLcmze] 



et al. (2008a I, the "cool clumps" would have to occupy 
70-90% of the space within 250 kpc radius, assuming that 
the two temperature phases are in pressure equilibrium. 
In conclusion , we find the scenario proposed by |Lemze| 
e t al.| ( |2008a| ) unhkely. 

ii'urther studying the ratio of Fe xxvi Ka and Fe xxv 
Ka emission lines, we conclude that these show no signs 
of a multi-temperature projection and the best fit of this 
ratio implies a single temperature consistent with the 
continuum temperatures from both XMM-Newton MOS 
and the Chandra data when the absorption edge is mod- 
eled. 

The discrepancy between lensing and X-ray mass es- 
timates remains, particularly in the r < 200 /i~^kpc re- 
gion. Our X-ray mass profile shows consis tent results 
compared to those from weak lensing (e.g., Broadhursi 
|et al. 2005bl [Lim ousin eF aljlSOOTl lUmetsu &: Broad 
nurst 2008; Corlcss et al. 2009). Strong lensing mass 



results (e.g., Broadhurst et al.|2005a 


Halkola et al. 12006 


Limousin et al.|2007 


) , but none ot them agrees with those 



Wu 2002 An 



dersson fc M adejski"'2004': 'Kiemer-S0 rensen et al.potigf :^ 
Using a simple ellipsoidal modeling of the cluster with 



simple ellips 

the major axis along the line of sight, we find that the 
projected mass, as derived from the X-ray analysis, in- 
creases by a factor of 1.6 assuming an axis-ratio of 0.6. 
We conclude that the mass discrepancy between lensing 
and X-ray derived masses can be alleviated by line of 
sight ellipticity and that this also can explain the high 
concentration parameter found in this cluster. 



We thank the anonymous referee for valuable sugges- 
tions on the manuscript. Support for this work was pro- 
vided by NASA through SAO Award Number 2834-MIT- 
SAO-4018 issued by the Chandra X-Ray Observatory 
Center, which is operated by the Smithsonian Astrophys- 
ical Observatory for and on behalf of NASA under con- 
tract NAS8-03060. EP sincerely thanks John Arabadjis 
for his support during early phases of this project. 



REFERENCES 



Allen, S. W., Rapetti, D. A., Schmidt, R. W., Ebeling, H., 
Morris, R. G., & Fabian, A. C. 2008, MNRAS, 383, 879 

Allen, S. W., Schmidt, R. W., & Fabian, A. C. 2001, MNRAS, 
328, L37 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 
197 

Andersson, K., Peterson, J. R., Madejski, G., & Goobar, A. 2009, 
ArXiv e-prints 

Andersson, K. E., & Madejski, G. M. 2004, ApJ, 607, 190 
Arabadjis, J. S., Bautz, M. W., & Arabadjis, G. 2004, ApJ, 617, 
303 

Arnaud, K. A. 1996, in Astronomical Society of the Pacific 
Conference Series, Vol. 101, Astronomical Data Analysis 
Software and Systems V, ed. G. H. Jacoby & J. Barnes, 17 — h 



Broadhurst, T., Bemtez, N., Coe, D., Sharon, K., Zekser, K., 
White, R., Ford, H., Bouwens, R., Blakeslee, J., Clampin, M., 
Cross, N., Franx, M., Frye, B., Hartig, G., lUingworth, G., 
Infante, L., Menanteau, F., Meurer, G., Postman, M., Ardila, 
D. R., Bartko, F., Brown, R. A., Burrows, C. J., Cheng, E. S., 
Feldman, P. D., Golimowski, D. A., Goto, T., Gronwall, C, 
Herranz, D., Holden, B., Homeier, N., Krist, J. E., Lesser, 
M. P., Martel, A. R., Miley, G. K., Rosati, P., Sirianni, M., 
Sparks, W. B., Steindling, S., Tran, H. D., Tsvetanov, Z. I., & 
Zheng, W. 2005a, ApJ, 621, 53 

Broadhurst, T., Takada, M., Umetsu, K., Kong, X., Arimoto, N., 
Chiba, M., & Futamase, T. 2005b, ApJ, 619, L143 

Cavagnolo, K. W., Donahue, M., Voit, G. M., & Sun, M. 2008, 
ApJ, 682, 821 

Cavaliere, A., & Fusco-Femiano, R. 1978, A&A, 70, 677 
Comerford, J. M., & Natarajan, P. 2007, MNRAS, 379, 190 



Discrepant Mass in A1689 



17 



Corless, V. L., King, L. J., & Clowe, D. 2009, MNRAS, 393, 1235 
Czoske, O. 2004, in lAU Colloq. 195; Outskirts of Galaxy 

Clusters: Intense Life in the Suburbs, ed. A. Diaferio, 183—187 
De Luca, A., & Molendi, S. 2004, A&A, 419, 837 
de Plaa, J., Werner, N., Bleeker, J. A. M., Vink, J., Kaastra, 

J. S., & Mendez, M. 2007, A&A, 465, 345 
Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 
Diego, J. M., Sandvik, H. B., Protopapas, P., Tegmark, M., 

Bem'tez, N., & Broadhurst, T. 2005, MNRAS, 362, 1247 
Faltenbacher, A., Kravtsov, A. V., Nagai, D., & Gottlober, S. 

2005, MNRAS, 358, 139 

Forman, W., Jones, C, Churazov, E., Markevitch, M., Nulsen, P., 
Vikhlinin, A., Begelman, M., Bohringer, H., Eilek, J., Heinz, S., 
Kraft, R., Owen, F., & Pahre, M. 2007, ApJ, 665, 1057 

Frye, B. L., Coe, D., Bowen, D. V., Bem'tez, N., Broadhurst, T., 
Guhathakurta, P., lUingworth, G., Menanteau, F., Sharon, K., 
Lupton, R., Meylan, G., Zekser, K., Meurer, G., & Hurley, M. 
2007, ApJ, 665, 921 

Gavazzi, R. 2005, A&A, 443, 793 

Halkola, A., Seitz, S., & Pannella, M. 2006, MNRAS, 372, 1425 
Hashimoto, Y., Bohringer, H., Henry, J. P., Hasinger, G., & 

Szokoly, G. 2007, A&A, 467, 485 
Hennawi, J. F., Dalai, N., Bode, P., & Ostriker, J. P. 2007, ApJ, 

654, 714 

Iwamoto, K., Brachwitz, F., Nomoto, K., Kishimoto, N., Umeda, 
H., Hix, W. R., & Thielemann, F.-K. 1999, ApJS, 125, 439 

Jeltema, T. E., Hallman, E. J., Burns, J. O., & Motl, P. M. 2008, 
ApJ, 681, 167 

Kaastra, J. S., & Mewe, R. 2000, in Atomic Data Needs for X-ray 
Astronomy, ed. M. A. Bautista, T. R. Kallman, & A. K. 
Pradhan, 161-+ 

Kaastra, J. S., Tamura, T., Peterson, J. R., Bleeker, J. A. M., 
Ferrigno, C., Kahn, S. M., Paerels, F. B. S., Piflaretti, R., 
Branduardi-Raymont, G., & Bohringer, H. 2004, A&A, 413, 415 

Kawahara, H., Suto, Y., Kitayama, T., Sasaki, S., Shimizu, M., 
Rasia, E., & Dolag, K. 2007, ApJ, 659, 257 

Kotov, O., & Vikhlinin, A. 2005, ApJ, 633, 781 

Kravtsov, A. V., Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128 

Lau, E. T., Kravtsov, A. V., & Nagai, D. 2009, ArXiv e-prints 

Leccardi, A., & Molendi, S. 2007, A&A, 472, 21 

— . 2008, A&A, 486, 359 

Lemze, D., Barkana, R., Broadhurst, T. J., & Rephaeli, Y. 2008a, 
MNRAS, 386, 1092 

Lemze, D., Broadhurst, T., Rephaeli, Y., Barkana, R., & Umetsu, 
K. 2008b, ArXiv e-prints 

Limousin, M., Richard, J., Julio, E., Kneib, J.-P., Fort, B., 
Soucail, G., Eh'asdottir, A., Natarajan, P., Ellis, R. S., Small, 
L, Czoske, O., Smith, G. P., Hudelot, P., Bardeau, S., Ebeling, 
H., Egami, E., & Knudsen, K. K. 2007, ApJ, 668, 643 

Lokas, E. L., Prada, F., Wojtak, R., Moles, M., & Gottlober, S. 

2006, MNRAS, 366, L26 

Markevitch, M., Bautz, M. W., Biller, B., Butt, Y., Edgar, R., 
Gaetz, T., Garmire, G., Grant, C. E., Green, P., Juda, M., 
Plucinsky, P. P., Schwartz, D., Smith, R., Vikhlinin, A., Virani, 
S., Wargelin, B. J., & Wolk, S. 2003, ApJ, 583, 70 



Mathiesen, B. F., & Evrard, A. E. 2001, ApJ, 546, 100 
Matsushita, K., Belsole, E., Finoguenov, A., & Bohringer, H. 

2002, A&A, 386, 77 
Mazzotta, P., Rasia, E., Moscardini, L., & Tormen, G. 2004, 

MNRAS, 354, 10 
Miralda-Escude, J., & Babul, A. 1995, ApJ, 449, 18 
Nagai, D., Vikhhnin, A., & Kravtsov, A. V. 2007, ApJ, 655, 98 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 

493 

Neumann, D. M. 2005, A&A, 439, 465 

Nevalainen, J., Markevitch, M., & Lumb, D. 2005, ApJ, 629, 172 
Oguri, M., & Blandford, R. D. 2009, MNRAS, 392, 930 
Oguri, M., Takada, M., Umetsu, K., & Broadhurst, T. 2005, ApJ, 
632, 841 

Piffaretti, R., Jetzer, P., & Schindler, S. 2003, A&A, 398, 41 
Pointecouteau, E., Arnaud, M., Kaastra, J., & de Plaa, J. 2004, 

A&A, 423, 33 
Pratt, G. W., & Arnaud, M. 2002, A&A, 394, 375 
Protassov, R., van Dyk, D. A., Connors, A., Kashyap, V. L., & 

Siemiginowska, A. 2002, ApJ, 571, 545 
Rasia, E., Ettori, S., Moscardini, L., Mazzotta, P., Borgani, S., 

Dolag, K., Tormen, G., Cheng, L. M., & Diaferio, A. 2006, 

MNRAS, 369, 2013 
Riemer-S0rensen, S., Paraficz, D., Ferreira, D. D. M., Pedersen, 

K., Limousin, M., & Dahle, H. 2009, ApJ, 693, 1570 
Saha, P., Williams, L. L. R., & Ferreras, I. 2007, ApJ, 663, 29 
Sanders, J. S., & Fabian, A. C. 2002, MNRAS, 331, 273 
— . 2007, MNRAS, 381, 1381 

Sarazin, C. L. 1988, X-ray emission from clusters of galaxies 

(Cambridge Astrophysics Series, Cambridge: Cambridge 

University Press, 1988) 
Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, 

J. C. 2001, ApJ, 556, L91 
Snowden, S. L., Mushotzky, R. F., Kuntz, K. D., & Davis, D. S. 

2008, A&A, 478, 615 
Tyson, J. A., Wenk, R. A., & Valdes, F. 1990, ApJ, 349, LI 
Umetsu, K., & Broadhurst, T. 2008, ApJ, 684, 177 
Vikhlinin, A. 2006, ApJ, 640, 710 

Vikhlinin, A., Burenin, R. A., Ebeling, H., Forman, W. R., 
Hornstrup, A., Jones, C, Kravtsov, A. V., Murray, S. S., 
Nagai, D., Quintana, H., & Voevodkin, A. 2009, ApJ, 692, 1033 

Vikhlinin, A., Forman, W., & Jones, C. 1999, ApJ, 525, 47 

Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C, Markevitch, 
M., Murray, S. S., & Van Speybroeck, L. 2006, ApJ, 640, 691 

Vikhlinin, A., Markevitch, M., Murray, S. S., Jones, C, Forman, 
W., & Van Speybroeck, L. 2005, ApJ, 628, 655 

Werner, N., Durret, F., Ohashi, T., Schindler, S., & Wiersma, 
R. P. C. 2008, Space Science Reviews, 134, 337 

Xue, S.-J., & Wu, X.-P. 2002, ApJ, 576, 152 

Zekser, K. C., White, R. L., Broadhurst, T. J., Bem'tez, N., Ford, 
H. C, lUingworth, G. D., Blakeslee, J. P., Postman, M., Jee, 
M. J., & Coe, D. A. 2006, ApJ, 640, 639 



APPENDIX 



INSTRUMENTAL ABSORPTION IN CHANDRA DATA 



To see whether the absorption feature found in A1689 (?2.4| is a cahbration problem, we analyzed other archived 
Chandra ACIS-I observations. Objects that have high-quality data and are relatively easy to model are pulsar wind 
nebula G021.5-00.9, elliptical galaxy NGC 4486 (M87), and the Coma cluster. Unfortunately, due to the high galactic 
absorption in G021.5-00.9, it is not suitable to use those observations to verify the instrumental absorption. Details of 



the datasets we used are listed in Table 14 These observations have low Nh (< 3 x 10 cm ) , low background level (< 



2% in 1.7-2.0 keV band), an d high signaPto-nois e ratios. Although the central part of M87 has very complex structures 



produced by the AGN (e.g., Forman et al.|2007[) , the XMM-Newton observation indicates that the intracluster medium 

;side those regions (iMatsushita et al. 20021. We extracted the spectra from 



is likely to be single-phase in nature outside 

r = 6' — 7.5' of M87 and r < 5' of Coma cluster for each ACIS-1 chip and ht with an absorbed single -temperature APEC 
mode l, multiplied by an absorption edge. The column density was fixed at the Galactic value (Dickey & Lockman 
1990 1. The redshift and all the elemental abundances, except Al, were free to vary. Parameters for this absorption edge 
are listed in Table[T5l For M87 whose emission is dominated by lines, these paramet ers are sensitive to the choice of 



the plasma model. Results from the latest MEKAL-based model, SPEX (version 2.0; Kaastra fc Mewe]|2000 1 , are also 
shown in Table 15 In general, an absorption at ~ 1.75 keV with an optical depth oi 0.1-0.15 is seen in the datasets. 
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TABLE 14 
Chandra Observation Log 



Name 


(102°cm-2) 


z 


ObsID 


Data 
Mode 


Obs. Date 


Exp. 

(ks) 


Background Norm. 
10 11 12 13 


Region 


(%) 


S/N'= 


M87 


2.59 


0.00423 


5826 


VFAINT 


2005-03-03 


125.5 


1.03 


1.00 


1.04 


1.07 


6'-7.5' 


1.5 


235 








5827 


VFAINT 


2005-05-05 


154.4 


1.10 


1.09 


1.09 


1.13 




1.5 


262 








6186 


VFAINT 


2005-01-31 


50.7 


1.02 


1.04 


1.04 


1.06 




1.5 


150 








7212 


VFAINT 


2005-11-14 


64.5 


1.15 


1.15 


1.13 


1.23 




1.7 


168 


Coma 


0.89 


0.0231 


9714 


VFAINT 


2008-03-20 


29.6 


1.32 


1.34 


1.18 


1.43 


<5' 


1.9 


134 



°- [Dickey fc Lockman| | |1990[ l 
^ Background fraction in 1.7-2.0 keV band. 
Signal-to-noise ratio in 1.7-2.0 keV band. 



TABLE 15 
Absorption edge parameters 



APEC 



SPEX 



Name ObsID ccd 



^thresh 

(keV) 



^thresh 

(keV) 



M87 


5826 





1 


47+0.02 
^' -0.02 





00+0.02 

^*-0.02 


(4.8±[i;^ 


X 10- 


3 


1 


47+0.02 
^' -0.02 





in+o.02 

-^-"-0.02 


(4 


0+0. 7^ 
"-0.7^ 


X 


10 


-3 






1 


1 


70+0.02 
'•^-0.03 





07+0.02 
-0.02 




X 10- 


3 


1 


7n+0.02 
'"-0.03 





00+0.01 

U«-0.02 


(4 


4+0. 4-) 


X 


10 


-3 






2 


1 


7c:+0.01 
'■^-0.02 





1 9+0.03 
-^-^-0.02 




) X 10 


-3 


1 


70+0.01 
'•-'-0.02 





1 9+0.02 
-^-^-0.02 


(8 


4+1. 2-, 


X 


10 


-3 






3 


1 


7Q+0.02 
'^-0.02 





09+0.02 
"^-0.02 


(5.3i«:^ 


X 


10- 


3 


1 


7c;+0.03 
'•-'-0.02 





00+0.02 
US-0.02 


(4 


^-0.3^ 


X 


10 


-3 




5827 





1 


7c:+0.02 
'■^-0.02 





1 1 +0.03 
^^-0.02 


(4.91^7 


X 


10- 


3 


1 


79+0.02 
'^-0.01 





11+0.01 
-^--^--0.02 


(3 


+ 1. 0-, 

''-o.s) 


X 


10 


-3 






1 


1 


7.1+0.02 
'•-'-0.02 





in+0.02 
^^^-0.02 


(5.5l?j 


X 


10- 


3 


1 


79+0.02 
'^-0.02 





09+0.01 
""-0.02 


(3 


1+0. 3-, 
-^--0.7' 


X 


10 


-3 






2 


1 


01 +0.02 
"-'--0.02 





1 1 +0.03 
^^-0.02 


(4.ll«:| 


X 


10- 


3 


1 


74+0.03 
'^-0.04 





00+0.02 
U*-0.02 


(2 


4+0.3-1 


X 


10 


-3 






3 


2 


-14+0.03 
-^-^-0.04 





Qg+0.02 
'^'-'-0.02 


(3.8l«j 


X 


10- 


3 


2 


1 0+0. 04 
^•^-0.06 





Qg+0.02 
""-0.02 


(2 


q+0.5-, 
''-0.3-' 


X 


10 


-3 




6186 





1 


70+0.06 
' ^-0.06 





09+0.03 




X 


10- 


3 


1 


7n+0.04 
'"-0.04 





09+0.02 
""-0.03 


(4 


0+0. 7n 
"-0.7^ 


X 


10 


-3 






1 


1 


70+0.03 
'•-'-0.04 





1 1 +0.04 
^^-0.03 




X 


10- 


3 


1 


71 +0.03 
'^-1.71 





1 1 +0.03 
^^-0.03 


(4 


4+0.8-1 


X 


10 


-3 






2 


1 


79+0.07 
'^-0.06 





1 1 +0.05 
^^-0.03 




X 


10" 


3 


1 


gg+0.03 
'^"-0.10 





10+0.02 
^"-0.02 


(7 


9+0. 4-, 
"-0.7/ 


X 


10 


-3 






3 


1 


74+0.03 
'^-0.03 





1 9+0.03 
^^-0.03 




X 


10- 


3 


1 


70+0.04 
'"-0.03 





1 1 +0.03 
^^-0.03 


(7 


0+0.3-) 
°-0.3^ 


X 


10 


-3 




7212 





1 


77+0.03 
' ' -0.02 





-1 -1 +0.02 
-^--^--0.04 




X 


10" 


3 


1 


(=9+0. 09 
"^-1.62 





00+0.02 
'J«-0.03 


(1 


0+3. 4n 
"-1.0/ 


X 


10 


-3 






1 


1 


01+0.21 

^-'--o.os 





in+0.04 
^"-0.04 


(5.8l«:^ 


X 


10- 


3 


2 


QQ+0.04 

^^"-0.04 





-in+0.03 
^"-0.05 


(2 


0+0. 4-, 
"-0.4^ 


X 


10 


-3 






2 




1-74/ 





m +0.03 


(6.41^7 


X 


10- 


3 




1.74; 





09+0.03 
"^-0.02 


(1 


9+1. 0-, 
^-1.0' 


X 


10 


-3 


Coma 


9714 





'-■'"-0.02 
1 7C.+0.05 


0.14t;!;* 


(4.6liJ;« 


X 10" 


2 


1 7g+0.03 

^- "^-0.02 


O-^^tVot 




X 10 


-2 






1 




/o 1+0.5 

'^•^■^-o.o 


X 10" 


2 


1 7C-+0.05 






X 10 


-2 






2 


, 7c:+0.02 
'-■'^-0.02 


11+0.04 
'-'•^^-0.02 


,0 9+0.7 
^•-'■^-0.3 


X 10" 


2 


1 7C-+0.03 
^- '•'-0.03 


10+°°5 
"■^"-0.02 


/o 9+0. 6n 


X 10 


-2 






3 


1 77+0.02 
' ' -0.02 


n -10+0.03 

U-i«-0.03 




X 10" 


2 


1 77+0.02 
'■ ' ' -0.02 
1 70+0.05 

-'-^ '•^-o.og 


Q 17 + 0.03 

"■^'-0.03 


/o 7+0.3-1 
V'^-'-0.5^ 


X 10 


-2 






3 


, 7C+0.06 
-^■'3-0. 06 


'-'•'^•^-0.02 


(4.8l«:^ 


X 10" 


3 


03+° °2 
"■"•^-0.01 


^9 0+0. 3n 
^^■•-'-0.8/ 


X 10 


-3 


A1689 


540 


13 


1 7r, + 0.05 
'^-0.03 


'-'•^'^-0.04 


n 190+0002 
U.iMU_Q (JQ2 




1 79+0.05 
'^-0.04 


"■""-0.03 




0.189: 


-0.004 
-0.003 






1663 


3 


, 7C-+0.04 
^■'■^-0.03 


'-'•^"-0.03 


(1 , 01 +0.006 

u.iyi_Q Qj^i 




1 7C-+0.04 
' •'-0.03 


16+°°'' 
"■J^"-0.04 




109+0.005 

"•1"^-0.012 






5004 


3 


^■'"-0.03 


15+0.02 
'-'•^■'-0.03 


n 109+0.002 




1 76+f-03 
J^- '0-0. 03 


15+0.03 
"■ ^^•'-0.03 




100+0.003 
0.1«U_o 






6930 


3 


1 74+0.01 
-'-■'*-0.02 


'-'•-'-^-0.01 


n ias+0.003 




1 7C-+0.01 
'•'-0.01 


"■^^-0.01 




100+0.002 
0-1OO-0.002 






7289 


3 


1 70+0.02 
'-■'■^-0.02 


-14+0.01 
'-'•-'-^-0.02 


n -loq+0.002 

U.iOO_Q (JQ3 




1 70+0.02 
'•^-0.02 


I'^+o.oi 

"■^^•^-0.01 




0.18lj 


-0.003 
-0.003 





However, for ObsID 7212, the absorption depth is determined less than 0.05 (Icr) on ACIS-I2 and 13. 



